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■Abstract 


Experimental  techniques  for  the  measurement  of  natural  extremely 
low  frequency  (EIF)  electromagnetic  noise  are  described,  including  methods 
for  absolute  calibration  of  magnetic  and  electric  field  sensors  in  this 
frequency  range  (3  to  30  Hz).  After  a  review  of  the  theory  applicable  to 
earth-ionosphere  cavity  resonances  a  method  is  developed  which  permits 
location  of  two  simultaneously  active,  major  thunderstoi®  regions  on  the 
surface  of  the  earth  by  analysis  of  Elf*  spectra.  In  particular  use  is 
made  of  electric  to  magnetic  field  ratios  and  of  ratios  involving  the 
power  spectnm  magnitudes  at  adjacent  resonant  peaks.  limitations  of  the 
method  are  analyzed,  particularly  those  related  to  f-®  background  noise 
(f  =  frequency)  and  to  quasi-uniformly  distributed,  minor  lightning  ac¬ 
tivity  which  adds  to  the  noise  originating  in  the  major  thunderstorm 
regions.  Appropriate  methods  of  spectral  analysis-appropriate  fre¬ 
quency  resolutions  and  integration  times— are  determined  and  the  method 
is  illustrated  by  experimental  data.  It  is  also  indicated  how  the  source 
location  procedure— in  conjunction  with  exact  data  (+_  0.25  Hz)  describing 
the  position  of  the  resonant  peaks— may  be  used  to  evaluate  Ionospheric 
conductivity  profiles  extending  downward  to  about  40  ka. 
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I.  INTRODUCTION  -  OBJECTS  AND  OUTLINE  OP  THE  STUDY  AND 


SUMMARY  OP  EXPERIMENTAL  METHODS 


1.1  Objects 

This  project  was  concerned  with  the  observation  and  explanation  of 
electromagnetic  noise  at  extremely  low  frequencies — approximately  Hz 

to  25  Hz — and  particularly  with  noise  near  the  Schumann  earth-ionosphere 
cavity  resonances  at  about  8,  14  and  20  Hz.  Objects  of  the  investigation 
were  the  identification  of  noise  sources,  the  identification  of  pertinent 
propagation  phenomena  and  the  identification  of  other  processes  which  de¬ 
termine  properties  of  ELF  noise  such  as  diurnal  variation  of  its  ampli¬ 
tude  and  of  its  frequency  spectrum. 

Possible  applications  of  the  observation  and  analysis  of  EIF  noise 
which  have  become  apparent  in  the  course  of  the  work  are  (l)  the  possi¬ 
bility  of  locating  the  position  of  major  thunderstorm  regions  on  the 
surface  of  the  earth,  (2)  identification  of  changes  in  the  effective 
conductivity  profile  (i.e.  electron  and  ion  density  and  collision  fre¬ 
quency  profile)  of  the  lower  Ionosphere  over  large  regions  and  (3)  pre¬ 
diction  f  some  types  of  Solar  Proton  events  and  associated  major  iono¬ 
spheric  disturbances. 

1.2  Outline  of  the  study 

The  major  part  of  the  present  report  is  devoted  to  a  detailed 
analysis  of  the  first  application — location  of  major  thunderstorm 
regions  (or  of  ether  ELF  sources)  by  analysis  of  ELF  observations 
at  a  single  station.  It  will  be  shown,  however,  that  this  location 
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procedure  is  also  fundamental  to  any  determination  of  ionospheric  con¬ 
ductivity  profiles,  since  accurate  information  about  such  profiles  cannot 
be  extracted  from  the  ELF  spectra  unless  the  relative  positions  of  source 
and  receiver  are  known  (Sections  V  and  X). 

Observations  indicating  substantial  deviations  from  normal  ELF  am¬ 
plitude  levels  and  normal  ELF  spectra  which  preceded  ten  Solar  Proton 
events  in  196?  and  1968  by  periods  between  eleven  rnd  twenty-three  hours 
will  not  be  covered  in  the  present  report,  because  they  are  not  well 
understood  and  because  the  number  of  observations  is  not  sufficient  to 
permit  definite  conclusions  without  further  observations.  Seme  possible 
source  and  propagation  mechanisms  were  analyzed  in  "Part  I"  of  this  final 
report  (APCRL-7O-OI43) . 

Since  much  of  the  information  obtainable  from  the  always  present 
natural  ELF  background  noise  is  contained  in  its  changing  power  spectrum, 
the  applicable  fundamentals  of  spectral  analysis  are  reviewed  in  Section 
II.  In  Section  III  the  characteristics  of  lightning  strokes  are  discussed, 
because  during  "normal"  periods — that  is  periods  when  ELF  noise  levels 
are  not  unusually  high  and  "Schumann"  resonances  are  clearly  noticeable— 
their  major  source  seems  to  be  lightning  activity. 

The  basic  theory  of  earth-ionosphere  cavity  resonances  due  to 
Schumann,  Wait,  Galejs,  D.  L.  Jones,  Madden  and  others  is  reviewed  in 
section  IV  while  details  of  a  new  ELF  source  location  procedure  are  dis¬ 
cussed  in  Sections  V  to  VII. 

The  effects  of  two  natural  phenomena  which  could  lead  to  large 
errors  in  the  source  location  procedure  are  analyzed  in  Section  VIII: 


(1)  Noise,  possibly  of  extra-terrestrial  or  Ionospheric  origin  and  having 

a  f”®  power  spectrum  (where  1  <  a  <  3)  is  always  present.  In  section  8.1 

it  is  indicated  how  large  this  background  noise  can  be  before  it  will 

prevent  extraction  of  source  information  from  ELF'  data.  (2)  The  source 

location  procedure  assunes  that  the  ma.lor  thunderstorm  activity  at  any 
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one  time  is  restricted  to  one  or  two  major  regions  such  as  a  10°  km  area 
(approx.  10  by  10  degrees)  between  1700  and  1900  local  time  at  an  equa¬ 
torial  land  mass.  Simultaneous  with  this  major  thunderstorm  activity 
there  are,  however,  many  smaller  thunderstorms  in  progress  over  the  rest 
of  the  earth’  surface.  Section  8.2  considers  the  effect  upo'  F  spectra 
(and  the  proposed  source  location  procedure)  of  various  possible  relative 
distributions  of  Dightning  activity  between  the  major  thunderstorm  centers 
and  the  more  uniformly  distributed  minor  lightning  activity. 

In  section  IX  the  source  location  procedure  is  illustrated  by  the 
analysis  of  experimental  ELF  power  spectra  and  in  section  X  conclusions 
are  drawn  and  directions  of  future  work  are  suggested.  In  particular  it 
is  indicated  how  the  source  location  procedure— employing  the  ratios  of 
the  power  spectrum  magnitudes  at  the  resonant  peaks — can  be  used,  in  con¬ 
junction  with  exact  data  (+_  0.25  Hz)  describing  the  position  of  the  reson¬ 
ant  peaks,  to  evaluate  Ionospheric  conductivity  profiles. 

1.3  Experimental  methods 

Observations  of  extremely  low  frequency  electromagnetic  noise  in  the 
earth-ionosphere  cavity  [Schumann]  resonance  region  are  made  on  the  Alton 
Jones  campus  of  the  University  of  Rhode  Island  [71°  44’  01"  W,  4l°  37* 
53"N]  in  West  Greenwich,  R.  I.  The  observation  site  is  in  a  rural  forest 


location,  far  from  sources  of  mannoade  electromagnetic  and  vibrational 
noise  (approximately  three  miles  from  the  nearest  public  road).  Varia¬ 
tions  of  the  horizontal  magnetic  field  and  of  the  vertical  electric  field 
in  the  frequency  region  between  approximately  3  and  40  Hz  (3  dB  points  of 
input  filter)  are  monitored  continuously;  the  vertical  electric  field 
data  cannot  be  used,  however,  during  heavy  local  precipitation.  Local 
lightning  activity,  when  it  occurs,  contaminates  both  electric  and  mag¬ 
netic  field  data,  but  thunderstorms  are  common  in  this  area  only  between 
May  and  September  and  the  annual  average  number  of  thunderstorm  days  is 
twenty-two. 

The  sensors  for  the  horizontal  magnetic  field  are  two  orthogonal  2 
meter  diameter  colls,  each  consisting  of  four  separate  sections  with  ap¬ 
proximately  11000  turns  of  number  30  wire.  Both  coils  are  electrostatically 
shielded  by  aluminum  housings  and  are  buried  in  the  ground  within  a  small 
mount  (approximately  2  meter  high  by  10  meter  in  diameter)  designed  so  as 
to  minimize  the  propagation  of  elastic  (acoustic)  surface  waves.  The  to¬ 
tal  output  from  each  coil  due  to  thermal  noise  within  the  coil  and  due  to 
pick-up  of  ambient  noise  is  approximately  4.5  microvolt  measured  through 
a  7.74  to  8.30  Hz  (one-tenth  octave)  filter  during  a  "quiet”  period  (no 
disturbances  on  the  "summary  record"  described  below  and  illustrated  by 
figure  1-5). 

The  magnetic  field  equipment  is  calibrated  at  least  once  each  week 
by  the  transmission  of  an  8  Hz  signal  from  a  calibrating  coil  for  a 
period  of  10  to  15  minutes.  The  calibrating  coil  produces  a  horizontal 
magnetic  field  intensity  of  exactly  10”^  ampere  (12.56  milligamma) 
at  the  location  of  the  receiving  coil  which  in  torn  produces  a  90  micro- 
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volt  output  voltage  in  the  7.74  to  8.30  He  band. 

The  output  of  each  receiving  coil  is  applied  to  an  amplifying  and 
filter  system  which  compensates  for  tlw  rising  anplltude-frequency 
characteristic  of  the  coil  so  that  the  output  of  the  system  is  propor¬ 
tional  to  the  magnetic  field  and  not  to  its  derivative.  The  frequency 
response  of  the  entire  receiving  system  [including  the  receiving  coll] 
is  flat  between  5  Hz  and  30  Hz  (+_  0.5  dBj  3  dB  points  are  at  3  Hz  and 
40  Hz)  for  constant  E- fie Id  input  (constant  current  in  the  calibrating 
coll).  The  output  of  the  simplifying  system  as  well  as  timing  signals  are 
recorded  AM  on  half -inch  magnetic  tape  running  at  a  speed  of  0.0375 
inches  per  second. 

The  vortical  electric  field  was  detected  by  a  vertical  rod  antenna 
until  August  1969  when  it  was  supplemented  by  a  spherical  antenna.  The 
latter  is  mounted  approximately  30  feet  (9.15  meter)  above  ground  level 
and  consists  of  a  silver  conductive  coating  painted  on  the  interior  of 
a  plexiglass  sphere.  Inside  the  sphere  is  a  high  impedance  preamplifier 
connected  to  the  inner  conductive  walls.  The  entire  assembly,  although 
mounted  on  top  of  a  mast,  cam  be  lowered  easily  to  the  ground  for  service 
and  repair. 

The  design  of  the  amtenna  is  based  upon  e  air  lie  r  work  by  Ogawa  (1966) . 
Mounting  of  a  metal  antenna,  which  is  electrically  isolated  from  ground, 
at  a  height  of  10  meters  largely  eliminates  spurious  signals  due  to  ion 
motion  which  is  particularly  pronounced  within  the  first  two  meters  above 
ground.  Furthermore  the  capacitance  to  ground  of  the  small  metal  sphere 
(radius  a  =  7»8l3  inches)  is  rather  insensitive  to  small  changes  in  height 
since  for  (a/h)  <  (l/lO)  its  value  is  given  by  C  »  2ncQa(2  +  a/h). 


Consequently  C  and  the  received  field  intensity  are  rather  insensitive 
to  lateral  motion  or  bending  of  the  supporting  mast  during  heavy  winds. 

The  antenna  differs  from  Qgawa's  "ball  antenna"  by  being  surrounded  with 
a  relatively  thick  (0.125  inch)  dielectric  coating  which  acts  as  a  "ra- 
doroe"  making  the  antenna  usable  even  during  moderate  precipitation. 

Since  absolute  calibration  of  the  electric  as  well  as  the  magnetic 
sensors  is  required  for  the  type  of  processing  described  later  in  this 
report  (for  example  the  effective  wave  impedance  |e|/|h|  is  of  interest), 
it  was  necessary  to  obtain  means  for  generating  a  known  standard,  ELF 
field  in  the  vicinity  of  the  receiving  antenna.  This  is  done  by  feeding 
a  2.86  meter  long  vertical  calibrating  antenna.  The  vertical  field  gener¬ 
ated  by  such  an  antenna  can  be  shown  to  be  given  by 


E_ 


where 


Vr  i_ 

4xw(H-f  \L 


(1.1) 


1^  =  base  current  of  the  transmitting  antenna 

H  =  elevation  of  the  top  of  the  transmitting  antenna  above  a  per¬ 
fectly  conducting  ground 

a  =  elevation  of  the  base  of  the  transmitting  antenna 

uj  =  radian  frequency 

z1  =  elevation  of  receiving  point 

rQ,  r^,  r^,  are  defined  by  figure  1.1 

This  expression  for  E  is  based  upon  the  equations  given  by  Jordan  (1950) 

z 

for  the  near-field  of  a  vertical  antenna.  In  deriving  (l.l)  appropriate 
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approximations  were  made  to  take  into  account  the  extremely  short  length 

of  the  antenna  and  the  extremely  short  observation  distance  in  eooparl- 

son  with  one  free -space  wavelength!  k  is  equal  to  3(10^)  meter  at  10  Be. 

It  is  apparent  from  (l.l)  that  for  an  accurate  determination  of  the 

ambient  field  in  the  vicinity  of  a  vertical  1 SIP  calibrating  antenna  the 

•  • 

only  quantity  which  must  be  measured,  in  addition  to  distances,  is  the 
base  current.  This  can  be  done,  Indirectly,  by  evaluating  the  capaci¬ 
tance  to  ground,  C  ,  of  the  transmitting  antenna.  The  antenna  current 

A 

is  then  equal  to  VuC  where  V  is  the  antenna  terminal  voltage. 

ft  ft  A 

Using  a  circuit  suggested  by  H.  L.  KSnlg  and  illustrated  in  Figure 
1.2  the  antenna  capacitance  is  given  by 


Ca  -  (R-l) (Cj  +  C2) 


where  R  is  defined  by 


R  = 


Z1  +  z2 


z 

P 


Za  Z2 

z  +  z 
a 


2 


(1.2) 


(1.3) 


(1.*) 


For  a  2.58  meter  calibration  antenna  used  on  the  West  Greenwich 

field  site  this  measurement  gave  C  =  22.2pf  while  a  calculation  based 

a 

upon  a  formula  given  by  Watt  (1967) 


C  « 
a 


log 


10 


24.l6(H-aj 

luprfl 


pf 


(1.5) 


rave  a  value  of  22.8  pf. 

In  (1.5)  the  quantity  (H-a)  is  the  antenna  length  as  shown  on  Fig¬ 
ure  1.2,  d  is  the  antenna  diameter  and  k  is  a  factor  depending  upon  the 


!*!*<* 


J/4':  .  ,  .  k..^.»  ,  _ 

i 
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ratio  of  antenna  elevation  above  ground  to  antenna  length.  For  the  value 
of  a  «  2.86  in  which  was  used  [a/(B-a)]  =  1.11  and  k  =  0.196. 

From  the  measured  capacitance  of  the  ■neceiviTfl  antenna  ,  which  is 
normally  used  for  the  reception  of  SIP  noise,  the  measured  amplifier 
gain  and  the  amplifier  output  voltage  one  can  evaluate  the  receiving  an¬ 
tenna  terminal  voltage.  If  the  receiving  antenna  terminal  voltage  thus 
calculated  agrees  with  the  terminal  voltage  obtained  by  integrating  over 
the  length  of  the  receiving  antenna  (for  a  linear  vertical  receiving  an¬ 
tenna)  the  field  due  to  the  calibrating  antenna  calculated  by  (1.1),  it 
becomes  possible  to  translate  receiving  antenna  voltage  accurately  into 
vertical  ELF  electric  field  values.  (For  a  spherical  receiving  antenna 
the  received  voltage  V  during  calibration  must  agree  with  the  calculated 
voltage  obtained  from  V  *  Q/C  where  Q,  is  found  by  integrating  the  calcu¬ 
lated  radial  component  of  the  calibrating  field  over  the  surface  of  thc- 

r__  _ 

receiving  antenna j  q  =  ejE°ds). 

t 

]  For  the  vertical  "whip"  ELF  receiving  antenna  which  has  been  used 

in  1967  and  1968  this  calibration  procedure  (at  10  Hz)  leads  to  the  fol¬ 
lowing  values  of  terminal  voltage: 

From  integration  of  the  radiated  calibrating  field  over  the  receiv¬ 
ing  antenna  -  23.96  mV 

From  amplifier  output  voltage,  amplifier  gain  and  receiving  antenna 
capacitance  -  25.89  mV 

It  appears  therefore  that  the  absolute  value  of  the  received  vertical 
electric  field  can  be  determined  with  a  precision  of  at  least  eight  per¬ 
cent.  Employing  the  described  calibration  procedure  one  finds  that  the 
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vertical  electric  field  at  the  Vest  Greenwich  field  site  hear  8  Hz  re¬ 
nal  ned  near  1.8  ailli-volt/neter/^S  on  a  typical  "quiet"  sunaer  day. 

Per  rapid,  visual  survey  of  field  intensity  variations  in  the  8  Be 
[first  Schumann  resonance]  region,  " stannary  records”  are  obtained  as 
follows*  the  nagnetlc  tape  is  replayed  in  the  laboratory  at  a  speed  of 
1.873  inches  per  second  [30  tines  the  recording  speed]  into  a  band-pass 
filter  having  3  dB  points  at  300  Hz  and  440  Hz,  corresponding  to  6  & 
and  8.8  Hz  in  the  original  recording.  The  output  of  this  filter  is 
passed  through  the  rectifier  and  smoothing  network  shown  in  Figure  1.3. 

The  response  of  this  filter  to  application  and  removal  of  a  constant 
400  Hz  signal  at  AA  is  shown  in  Figure  1.4.  The  charging  tine  constant 
Tq  of  the  filter  [time  needed  fo*'  output  to  reach  1  -  e-1  of  final  value] 
is  7-5  seconds.  The  measured  discharge  tine  constant  T^  is  31  seconds; 
the  discharge  function  actually  has  the  fora 

1  r  “t/^P-  — t/’i  p  -« 

v  -  [Tx  -  T2]_1  [Txe  1  -  T2e  2  ]  , 
but  since  »  T2  it  follows  that 
-tAi 

V  »  e  .  [T2  =  1.6  seconds] 

Since  the  input  time  constant,  Tq  =  7.3  seconds,  corresponds  to 
[7-5(50)/6o]  =  6.25  minutes  of  recording  time,  fluctuations  in  the  input 
signal  which  are  appreciably  shorter  than  six  minutes  will  not  appear  on 
the  paper  record  which  is  produced  by  the  output  of  the  network  in  Fig¬ 
ure  1.3.  Furthermore*  low  level  fluctuations  which  occur  within  a  period 
of  T  [50]/60  =  42  minutes  after  a  sustained  large  signal  was  received  will 
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not  appear  on  the  output  chart.  The  "avranary  records'*  illustrated  by 
Figure  1.5  (and  regularly  published  in  Section  5  of  the  AFCRL  "Geo¬ 
physics  and  Space  Data  Bulletin")  should  therefore  be  used  primarily  to 
locate  periods  of  sustained  unusual  activity)  for  such  periods  the 
original  wide -band  [jj  Hz  to  40  Hz]  magnetic  tape  recordings  may  then 
be  examined  in  more  detail.  The  "sxmnary  records"  can  also  be  used  for 
comparison  with  thunderstorm  data  and  records  of  solar  and  geomagnetic 
activity  to  establish  which  events  have  significant  effect  rqxm  magnetic 
field  variations  in  the  8  Hz  region. 

Calibration  signals  are  indicated  by  C  on  the  records;  periods  of 
man-made  interference  or  missing  data  are  indicated  by*-*;  local  thunder¬ 
storms  (within  about  25  miles  of  the  field  station)  are  shown  by  LS. 

_5 

The  10  amp  m  “  8  Hz  calibrating  signal  saturates  the  amplifier  of  the 
paper  chart  recorder  for  the  gain  setting  which  is  employed  (although  it 
is  10  dB  below  saturation  level  on  the  magnetic  tape) ,  it  can  therefore 
not  be  used  to  establish  relative  signal  levels,  except  when  records  are 
reproduced  at  half  gain.  Time  on  the  records  is  EST  with  heavy  vertical 
lines  indicating  0000  hours  EST. 

Fundamental  considerations  important  in  the  processing  for  detailed 
spectral  analysis  of  the  wide  band  magnetic  field  and  electric  field  data 
which  are  recorded  on  magnetic  tape,  are  discussed  in  Chapter  II  and  in 
Appendix  A-l.  Spectral  analysis  was  performed  by  either  digitizing  the 
data  and  employing  a  "Fast-Fourier-Transfom"  program  on  an  IBM  360/50, 
or  more  recently— achieving  a  reduction  of  analysis  time  by  a  factor  of 
about  JO/l— by  employing  a  commerical  ("3AIC0R")  special  purpose  analog- 
digital  system  for  Fourier  spectrum  analysis. 


H.  SPECTRAL  ANALYSIS 


The  convenience  of  the  frequency  concept  in  the  analysis  of 
physical  phenomena  is  well  known.  In  order  to  facilitate  the  compari¬ 
son  of  theoretical  with  experimental  results,  however,  a  method  of  in¬ 
vestigating  experimental  data  in  the  frequency  domain  is  needed.  This 
is  the  tool  of  spectral  analysis.  Determining  the  frequency  content  of 
a  time  signal  is  by  no  means  an  easy  task  to  accomplish  and  in  the  fol¬ 
lowing  paragraphs  some  of  the  problems  that  occur  in  spectral  analysis, 
as  well  as  methods  for  treating  them  are  discussed. 

The  frequency  content  of  a  time  function  is  given  by  its  Fourier 
transform.  A  time  signal,  x(t),  and  its  Fourier  transfom,  1(f)  are 
related  by  _ 

X(f)  =  j  x(t)e~i2ntt  dt  (2.1) 

JO 

and  +• 

x(t)  «  j  X(f)e127lft  df  (2.2) 

JO 

provided  certain  mathematical  conditions  on  x(t)  and  X(f)  are  satisfied 
[Weinberger,  1965].  These  conditions  are  always  satisfied  by  physically 
realizable  waveforms. 

Although  (2.1)  does  give  the  frequency  content  of  the  signal,  the 
result  is  complex,  i.e.  it  includes  the  phases  of  the  spectral  components. 
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Usually,  however,  phase  information  is  not  of  interest  and  a  positive, 
real  quantity  jl(f)}  is  used.  The  power  spectral  density,  Sx(f)  can 
be  defined  as 

S  (f)  *  liin  « 
x  T  -  ® 


I 


-T 


x(t)e“i2xftdt 


(2.5) 


1 

1 

| 

5 


5 

| 

I 


* 

; 
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It  is  called  a  power  density  because  Sx(f)  df  gives  a  measure  of  the 
power  in  the  frequency  band  between  f  and  f+df.  An  alternative  defini¬ 
tion  of  power  spectral  density  in  terns  of  the  autocovariance  function 
Hx(t)  is  sometimes  used.  The  autocovariance  function  is  defined  as 

T 

\(t)  -  11m  x(t)  x  (t+T)dt  (2.k) 

T  -  -  _T 

The  power  spectral  density  can  then  be  expressed  as  the  Fourier  trans¬ 
form  of  the  autocovariance  function. 

+  • 

Sx(f)  =  J  ^(Tje-211^  dT  (2.5) 

mt  ® 


The  definitions  (2.5)  and  (2.5)  are  completely  equivalent .  Es¬ 
timation  of  the  power  spectrun  by  (2.5),  which  until  recently  was  the 
most  commonly  used  method  for  estimating  the  power  spectrum,  involves 
two  computational  steps:  the  evaluation  of  the  autocovarlance  and  then 
the  evaluation  of  the  Fourier  transform  of  the  autocovariance.  In  the 
computation  of  the  autocovariance,  however,  lags  of  less  than  or  equal 
to  only  about  ten  percent  of  the  total  record  length  need  be  considered 
(Blackman  and  Tukey,  1958).  For  this  reason  the  evaluation  of  the  in¬ 
tegral  in  (2.5)  is  much  simpler  than  the  evaluation  of  the  integral  in 
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(2.3).  In  the  past,  significant  computer  time  savings  were  realized 
by  using  (2.5)  since  the  calculation  of  the  numerical  Fourier  transform 
was  a  lengthy  process  which  could  be  considerably  shortened  by  reducing 
the  domain  of  the  function  to  be  transformed. 

Cooley  and  Tukey  (1965)  have,  however,  developed  a  highly  efficient 
method  for  numerically  evaluating  a  Fourier  transform.  With  the  advent 
of  this  algorithm,  the  emphasis  has  shifted  to  using  (2.3)  for  the  evalu¬ 
ation  of  power  spectra.  In  fact,  the  Cooley-Tukey  algorithm  or  "Fast 
Fourier  Transform"  (FFT)  is  so  fast  that  even  if  the  autocovariance  is 
desired,  it  is  more  efficient  to  find  the  power  spectrum  by  the  FFT  and 
then  to  take  the  inverse  Fourier  transform  via  the  FFT,  yielding  the 
autocovariance.  The  power  spectral  density  can  be  efficiently  estimated 
by  evaluating  (2.3)  on  a  digital  computer  using  the  FFT. 

2.1  Resolution  and  Bias  of  the  Estimates 

Evaluation  of  equation  (2.1)  requires  knowledge  of  the  entire 

history  of  signal,  past  and  future.  This  is,  of  course,  not  possible 

and  we  must  use  an  approximation  of  the  form 

T/2 

X  (f)  =  f  x(t)e'l27lft  dt  (2.6) 

a  <i 

-T/2 

where  [-T/2,  T/2]  is  the  time  interval  of  the  signal  which  is  available. 
Equation  (2.6)  can  be  written  as 

+o= 

X  (f)  -  1  h(t)  x(t)e“i2nft  dt  (2.7) 

a 

—oc 

if  we  define  h(t)  by 


1;  -  T/2  £  t  £  T/2 


(2.8) 


« 


a(t)  =  { 

0;  t  <  -  T/2  or  t  >  T/2 

h(t)  is  often  called  a  "data  window"  since  it  determines  what  portion 
of  the  data  is  "seen"  by  our  analysis.  Using  (2.2),  we  can  express 
x( t)  in  terms  of  its  Fourier  transform  and  interchange  the  order  of 
integration  in  the  resulting  iterated  integral,  to  obtain 

•po 

Xa(f)  ,  J  X(f)  H(f-f’)  df*  (2.9) 

The  function  H(f),  called  a  "spectral  window",  is  the  Fourier  trans¬ 
form  of  h(t),  the  data  window,  and  is  shewn  in  Fig.  2.1.  The  effect 
of  employing  (2.6)  to  estimate  X(f)  is  clear  from  (2.9);  X  (f)  is  a 

a 

weighted  average  of  the  actual  Fourier  transform.  Since  the  width  of 
the  main  lobe  of  H(f)  is  proportional  to  T-1  (Fig.  2.1),  we  can  increase 
xhe  resolution  of  our  estimates  by  increasing  the  length  of  signal  that 
we  analyze.  In  fact,  as  T  becomes  infinite  (or  equivalently,  h(t)  =  1 
for  all  time)  H(f-f')  becomes  ^(f-f?)>  the  unit  impulse  function, 
(Papoulis,  1962)  and  X  (f)  =  X(f). 

u 

The  estimates  of  X(f)  can  be  biased  by  "leakage"  from  adjacent 
frequency  bands  due  to  the  side  lobes  of  H(f).  This  effect  becomes 
especially  troublesome  if  the  power  spectrum  of  x(t)  contains  strong 
sinusoidal  components  cr  large  peaks  (Hinich  and  Clay,  1968).  Bendat 
and  Piersall  (19 66)  have  approximated  this  additive  bias  of  spectral 
sstimates  as 
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where  Af  is  the  frequency  resolution  in  hertz  of  the  estimates.  That 
this  bias  term  is  negligible  for  earth-iono sphere  resonant  cavity  power 
spectra  is  shown  in  Appendix  A. 

2.2  Variability  of  the  Estimates 

To  numerically  estimate  the  Fourier  transform  of  a  time  function, 

the  integral  in  equation  (2.6)  must  be  calculated.  An  approximate 
* 

formulation  is 

N-l 

-i2imfnkAt 

An  =  A(nfQ)  =  At  )  x(kAt)c  1 

k=0  n  =  0,  1,  2,  ...N-l 

(2.10) 

where  At  is  the  tine  interval  between  the  samples  of  the  analog  to 
digital  converter  and  fQ  is  the  fundamental  frequency,  fQ  =  l/T  = 

1/NAt.  From  (2.3)  we  see  that 

Sn  =  T  KJ2  ;  n  «  °»  It  2.  ...N-l  (2.11) 


are  estimates  of  the  power  spectrum.  A  plot  of  versus  n  is  called 
the  raw  periodogram.  It  will  be  shown  that,  unless  modified,  the  raw 
periodogram  is  not  a  statistically  consistent  estimate  of  the  power 
spectrum. 

Suppose  that  a  frequency  resolution  of  A  f  hertz  is  desired  in  the 

power  spectrum.  To  achieve  this  resolution  M  +  1  values  of  S  could 

n 

be  averaged  as 


+M/2 

P  . . i—  r 

n  M  +  1  L, 


S. 

l+n 


(2.12) 


i=  -M/2 


*For  convenience,  we  have  shifted  the  time  interval  to  [o,  T], 


where  M  i r  use  nearest  wen.  integer  to  (Sf/f^).  The  resulting  would 
constitute  sa  estimate  of  the  power  spectrum  at  f  *  nf.  Hertz. 

V 

An  expression  for  the  variance  of  P_  will  he  derived  be  lev.  We 

n 

3?iall  aj8u,.A  that  the  tire  signal  is  a  member  of  a  white  stationary 
gauss ian  random  process,  or,  in  other  words,  that  'ach  tine  sanple 
x(kAt)  is  a  normally  distributed  random  variable  whose  statistics  are 
independent  of  k  and  that  the  power  spectrum  of  x(t)  is  a  constant  for 
all  frequencies.  Using  the  notation  a  x(kAt),  equation  (2.10)  can 
be  written 


A  *  6t 

n 


N-l 


k=0 


-  §SL  ku 
K 


(2.13) 


Equation  (2,13)  is  of  the  fora 


where 


b  a 

+  ib 

ft 

n 

R-l 

=  It 

T  \ 

k==0 

N-l 

br  -  -At  V  ^  sin  (  jp-  kn) 
k=0 


If  is  assumed  to  be  a  normally  distributed  random  variable,  the  an 
and  will  also  be  nonrally  distributed,  since  they  represent  linear 
combinations  of  normally  distributed  random  variables.  (Mood  and  Gray- 
bill,  1963).  It  will  also  be  assumed  that  E(x^)  =  0,  i.e.  the  signal  is 
passed  through  a  filter  or  some  device  that  removes  the  d-e  level.  In 


digital  spectral  analysis,  the  signal  is  usually  filtered  to  prevent 


aliasing  (Blackman  and  Tukey,  1938)*  If  the  gain  on  the  data  gathering 
equipment  is  adjusted  we  can  force 

Var(an)  -  Var(t>n)  •  1|  n  «  1,  2,  3,... N-l  (2.14) 

so  that  the  quantities 

|Aj2  =  an?  +  bn2;  n  =  1,  2,  ...N-l  (2.15) 

are  distributed  chi-square  with  two  degrees  of  freedom.  (Hood  and 
Graybill,  1963).  It  is  clear  that  the  mean  value,  (t  is  given  by 

H  =  E[jAn|2]}  a-  1.  2,  3t  ...N-l 

2 

and  that  the  variance,  a  is  given  by 

o2  =  E[(|Anj2  -  n)2];  n  =  1,  2,  3,  ...N-l 

2 

where  E  denotes  expected  value.  Since  a  is  a  measure  of  power,  the 

last  equation  is  a  result  of  our  white  noise  assumption. 

The  variance  of  a  chi-square  variate  with  two  degrees  of  freedom 

equals  its  mean.  (Mood  and  Graybill,  196 3) -  Thus, 

14/2 

E(Pn>  -  5iTT  I  E<KJ2)'“  <2‘16> 

i  «  -M/2 

The  random  variables  |A  |  with  n  =  1,  2,  3>  ...N-l  are  independent 
since  they  are  obtained  by  an  orthogonal  transformation  and  represent 
different  "components"  of  x(t).  The  variance  of  a  sum  of  independent 
random  variables  is  the  sum  of  the  variances  (Mood  and  Graybill,  1963). 


HI.  THE  LDHUKDO  STOCKS  AS  A  SOURCE  OP  SIP  WAVES 

3.1  The  Lightning  Discharge  Process 

Upon  application  of  an  electric  field,  the  constituent  atoss  sad 
molecules  of  a  gas  becone  polarised.  If  the  electric  field  is  increased 
until  the  ionisation  potential  is  achieved,  electrons  will  be  liberated 
from  their  parent  atone  and  be  accelerated  by  the  applied  field.  Under 
certain  conditions  (Loeb,  1939),  the  free  electrons  can  achieve  the 
momentum  necessary  to  ionise  other  atoms  by  collision.  The  fast  nerv¬ 
ing  electrons  can  also  excite  neutral  atoms  causing  them  to  emit  radia¬ 
tion  which  can  photoianize  other  atoms.  Hence,  electron  multiplication 
or  electron  avalanche  can  occur  in  an  ionized  gas.  The  lightning  dis¬ 
charge  is  an  exanple  of  this  phenomenon. 

The  lightning  flash  is  a  very  complex  mechanism  and  only  the  gen¬ 
eral  features.  Important  to  the  present  discussion,  will  be  presented. 

If  adjacent  regions  of  air  are  at  different  temperatures  convection 
takes  place:  warmer  air  tends  to  rise  and  cooler  air  descends.  In  a 
mature  thunderstorm  cell,  convection  is  occurring  on  a  large  scale 
and  a  strong  updraft  is  present.  It  is  this  tqxiraft  which  carries 
positive  charges  to  the  top  of  the  cloud.  Negative  charges,  which  are 
believed  to  reoide  on  heavier  particles  remain  in  the  lower  regions  of 
the  cloud.  The  negative  charges  in  the  bottom  of  the  thundercloud  in¬ 
duce  a  positive  surface  charge  on  the  conducting  earth  below.  The 


electric  field  in  the  region  between  the  lower  cloud  and  the  earth  in¬ 
creases  as  the  charge  separated  by  the  updraft  increases.  When  the 
breakdown  potential  is  reached,  the  air  below  the  cloud  becomes  ionised. 
The  ionization  will  first  occur  where  the  electric  field  ie  strongest, 
which  is  usually  just  below  the  bottom  of  the  cloud.  The  free  electrons 
in  the  ionized  air  are  then  accelerated  towards  the  earth  by  the  repul¬ 
sive  force  of  the  large  accumulation  of  negative  charge  in  the  lower 
portions  of  the  cloud;  any  motion  of  the  positive  and  negative  ions  can 
be  neglected  since  their  large  mass  renders  them  Immobile  relative  to 
the  electrons  (Schonland;  1964).  The  movement  of  the  electrons  takes 
place  in  a  series  of  steps  of  about  30  meters  in  length.  Hence,  this 
first  part  of  the  lightning  stroke  is  called  the  stepped  leader.  When 
the  stepped  leader  reaches  the  earth,  a  conducting  path  exists  between 
the  large  accunulations  of  charge  in  the  cloud  and  the  earth.  The  elec¬ 
trons  in  the  lower  part  of  the  leader  flow  into  the  ground  first,  then 
the  electrons  behind  them  so  the  process  appears  to  move  up  towards  the 
cloud.  The  return  stroke  is  actually  an  upward  moving  increase  in  the 
downward  velocity  cf  electrons.  It  is  the  return  stroke  that  produces 
the  intense  light,  heat,  arid  acoustical  energy  of  the  typical  lightning 
discharge. 

It  should  be  noted  that  the  lightning  discharge  to  earth  consists 
visually  of  several  leader-return  processes  in  succession,  with  the  av¬ 
erage  number  being  three;  up  to  42  strokes  in  a  single  discharge  to  earth 
have  however  been  observed  (Schonland,  1955)*  The  leaders  for  strokes 
subsequent  to  the  first  are  not  stepped,  but  move  continuously  towards 
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the  ground.  These  leaders  are  called  dart  leaders  to  distinguish  then 
from  the  slower  moving  stepped  leaders.  The  average  period  of  time  be¬ 
tween  strokes  is  about  0.05  secs. 

The  stepped  leader  travels  at  an  effective  speed  of  JQO  kn/sec. 

The  dart  leader  and  return  stroke  travel  at  speeds  of  roughly  2,000  and 
10,000  km/sec  respectively.  The  return  stroke  current  reaches  a  peak 
value  of  about  30,000  amps  within  6  usees  after  the  leader  reaches  the 
ground.  The  current  decays  to  half  the  peak  value  in  an  average  of 
24  usees  (Chalmers,  1967). 

Although  only  the  discharge  to  earth  has  been  discussed,  lightning 
flashes  that  never  reach  the  ground  are  quite  common.  These  discharges 
do  not  manifest  a  return  process  but  are  only  of  the  leader  type  (Schon- 
land,  1953)*  A  stroke  between  two  different  charge  concentrations  within 
the  same  cloud,  or  sheet  lightning  as  it  is  caanonly  called,  is  observed 

* 

frequently.  Air  discharges,  that  is,  lightning  flashes  passing  from  a 
thundercloud  into  clear  air  (Schonland,  1956)  as  well  as  strokes  from 
the  highest  regions  of  the  cloud  up  to  the  electrosphere  (workers  in 
atmospheric  electricity  refer  to  the  cosmic  ray  produced  ionized  layer 
as  the  electrosphere  and  to  the  higher  conducting  layers,  a  result  of 
solar  radiation,  as  the  ionosphere)  have  been  reported  (Chalmers,  1967)* 
Strokes  along  a  path  between  two  clouds  (intercloud  strokes)  are  also 
possible. 

The  characteristics  of  individual  lightning  discharges  probably 
are  determined  to  a  large  extent  by  many  climatological,  topographical, 
and  geographical  conditions.  For  example,  evidence  indicates  that  the 
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percentage  of  discharges  to  earth  depend*  upon  geographical  latitude 
(Chalmers,  1967).  Tropical  discharge  may  differ  significantly  from 
those  occurlng  in  temperate  zones  (Vatt ,  i960).  According  to  Schon- 
land  (1964) ,  "giant"  lightning  flashes,  with  pe»k  currents  greater  than 
200,000  amperes  have  only  been  observed  in  the  Swiss  Alps.  lightning 
strokes  to  the  Rnpire  State  Building  have  been  observed  to  begin  with  a 
leader  moving  up  towards  the  clouds. 

The  large  differences  in  the  properties  of  individual  lightning  dis¬ 
charges  constitute  a  possible  limitation  of  the  methods  to  be  reported 
here.  It  will  be  shown  below,  however,  that  many  of  these  difficulties 
can  be  overcoats  by  enploying  the  proper  characteristics  of  received  SI F 
signals  for  the  location  of  sources. 

3.2  The  lightning  Stroke  aa  a  Source  of  Radio  Waves 

Since  the  lightning  flash  consists  of  a  time  changing  flow  of  charged 
particles,  it  must  radiate  electromagnetic  energy.  In  view  of  the  large 
wavelength  of  ELF  radiation,  the  lightning  stroke  can  be  modeled  as  an 
infinitesimal  electric  dipole  of  current  dipole  moment  M(t)  eoul-meter/sec. 

According  to  Watt  (i960)  fields  radiated  by  leader  processes  in 
lightning  discharges  at  frequencies  below  5  kHz  are  not  significant.  One 
might  expect  this  to  be  true  for  the  stepped  leaders  which  are  relatively 
short  in  length,  but  this  may  not  be  the  case  for  the  dart  leaders  which 
are  as  long  as  the  return  stroke  and  may  have  currents  of  the  order  of 
10,000  amps  (Schonland,  1956).  Pierce  (1963)  has  formulated  an  empiri¬ 
cal  expression  for  the  time  rate  of  change  of  lightning  dipole  moment 


a  cloud  to  ground  stroke  as 


M(t)  -  80  h1exp(-h12t2)  +  ^^(-hgV) 

♦  35  h^ejp(-h^2t2) 


(5.1) 


where  the  tens  la  h^  aooouats  for  leader  prooeaaea,  the  tem  in  h^  ac- 
oouzxta  for  the  return  stroke,  and  h^  accounts  for  any  continuing  current 
after  the  return  stroke.  Por  a  cloud  die  charge,  only  the  first  tem 
applies.  The  spectral  composition  of  the  discharge  moment  can  be  easily 
shown  from  (3.1)  to  be 


M(w)  -  80  exp  (-^A^2)  ♦  5  exp  (-^Ahg2) 

2  o  (3.2) 

+  35  exp  (hu  A^j  )  coul-km 

The  values  of  1^,  h2  and  given  by  Pierce  (1963)  are  40,  40  x  10^,  and 
1(P.  Figure  (3.1),  which  shows  each  term  in  equation  (3.2)  as  well  as 
the  total,  M(w),  indicates  that  the  leader  is  much  more  significant  for 
the  first  three  Schumann  resonances  than  the  return  stroke.  It  is  clear 
that  there  is  disagreement  about  which  component  of  the  discharge  pro¬ 
cess  is  most  important  at  EIF. 

Lightning  strokes  can  have  both  vertical  and  horizontal  components. 
Several  workers  have  determined  from  theoretical  and  experimental  inves¬ 
tigations  that  it  is  the  vertical  currents  which  are  the  major  source  of 
low  frequency  (EIF  and  VLF)  electromagnetic  waves  in  the  earth  ionosphere 
cavity.  For  example  Williams  (1959)  has  shown  by  extensive  statistical 


studies  of  thunderstorm  data  that  the  major  sources  of  VLF  waves  are 
vertical  lightning  currents.  Scnmerfeld  (1949)  has  analytically  treated 
the  problem  of  radiation  from  both  vertical  and  horizontal  dipoles  over 
a  flat,  perfectly  conducting  earth  assuming  a  homogeneous  atmosphere. 

His  results  show  that  when  the  dipole  height  to  wavelength  ratio  is 
very  null  (as  is  the  case  for  ELF) ,  the  power  radiated  from  a  horizon¬ 
tal  dipole  is  many  orders  of  magnitude  less  than  the  power  radiated  by 
a  vertical  dipole  of  equal  moment.  Employing  a  more  realistic  model, 
Wait  (1962)  has  shown  that  the  ratio  of  the  vertical  electric  field  from 
a  horizontal  dipole  to  the  vertical  electric  field  from  a  vertical  di¬ 
pole  is  proportional  to  u c/a  which  for  o  =  10”^  mho/meter  is  equal 

D  © 

to  about  lO-'’  in  the  lower  Schumann  band.  The  problem  of  horizontal 

dipole  excitation  of  the  earth-ionosphere  waveguide  has  been  treated  in 

a  recent  paper  by  Galejs  (1968a).  The  horizontal  dipole  excites  both 

TM  and  TE  modes  but  all  modes  except  the  zero  order  TM  are  evanescent 

at  ELF.  The  vertical  dipole  is,  however,  a  much  more  efficient  ercitor 

of  the  TM  mode  than  the  horizontal  dipole, 
o 

The  reason  for  the  poor  radiation  properties  of  the  horizontal  di¬ 
pole  seems  to  be  intuitively  clear.  In  order  to  satisfy  the  boundary 
conditions  at  the  surface  of  the  earth,  the  '‘image”  of  the  source  must 
be  oppositely  directed.  If  the  source  and  its  image  are  very  close  (in 
units  of  wavelength  they  are  practically  coincident  at  ELF)  their  fields 
counteract  to  greatly  reduce  their  effective  radiation.  In  view  of 
these  considerations  and  the  experimental  evidence  it  is  generally  ac¬ 
cepted  (Harris  and  Tanner,  1962:  Galejs,  1961}  Raeraer,  1961)  that  only 
the  vertical  components  of  lightning  current  are  significant  in  the 


26 


production  of  terrestrial  EIFnoiao. 

3.3  Frequency  Spectrum  of  lightning  Stroke  Dipole  Moment 

There  are  very  few  data  available  on  ELF  lightning  spectra  (Klmpara, 
1963} .  Harris  and  Tanner  (1962)  have  estimated  the  frequency  variation 
of  a  quantity  "very  nearly  proportional  to  the  frequency  spectrum  of 
the  square  of  the  current  moment”  of  the  lightning  stroke  using  the  ELF 
data  of  Balser  and  Vagner  (i960).  Rycroft  (1963)  calculated  the  fre¬ 
quency  spectrum  of  lightning  from  plausible  time  variations  of  the  stroke 
dipole  moment  and  Raemer  (1961)  estimated  the  power  spectrum  of  a  median 
return  stroke  enqploying  the  lightning  stroke  statistics  compiled  by 
Williams  (1939).  Galejs  (1961)  approximated  the  Raemer  spectrum  by  the 
relation 

g(w)  o  const,  exp  [-9.1  x  10"3u]  (3*5) 

The  agreement  among  the  three  lightning  spectra  plotted  in  Fig, 

(3.2)  is  surprisingly  good  considering  that  completely  different  methods 


were  used  in  their  derivation, 


TV.  CAVITY  FIELDS  FROM  VERTICAL  ELECTRIC  DIPOLE  SOURCES 


4.1  The  Earth- Ionosphere  Cavity  Model 

The  conductivity  of  the  ionosphere  is  orders  of  magnitude  less 
than  the  effective  ground  conductivity  so  that,  in  the  ELF  range,  the 
ionosphere  is  the  major  contributor  to  wave  attenuation  and  the  earth 
can  be  assumed  perfectly  conducting  (Galejs,  1964a). 

Ionosphere  models  (Galejs,  1961;  Jones,  1967)  admitting  only  the 
vertical  inhomogeneity  lead  to  results  which  are  in  closer  agreement 
with  experiment  than  a  sharply  bounded,  cons cant  parameter  model. 
Horizontal  variations  in  ionospheric  parameters  are  usually  considered 
less  important  than  vertical  variations;  but  the  effects  of  lateral 
ionospheric  inhomogeneities  have  been  examined  by  Madder.-Thongjson  (1965) 
and  by  Nelson  (1967).  One  important  horizontal  variation  in  ionospheric 
properties  occurs  at  the  boundary  between  the  day  and  night  hemispheres; 
this  problem  has  been  treated  by  Galejs  (1964a),  but  his  method  is 
applicable  only  if  the  source  is  located  along  the  day-night  bcundary. 

The  effect  of  the  terrestrial  static  magnetic  field  on  the 
propagation  of  ELF  and  VLF  waves  has  recently  received  considerable 
attention  in  che  literature.  Galejs  (1968b)  has  shown  that  there  is 
very  little  change  in  the  attenuation  and  phase  velocity  of  JO  Hz  waves 
when  the  propagation  direction  is  east  to  west  or  west  to  east  for  a  45° 
magnetic  dj.p  angle.  For  other  propagation  directions  the  effect  of  the 


magnetic  field  is  even  smaller  (Galejs,  1968c).  The  values  of  attenua¬ 
tion  and  phase  velocity  computed  by  Galejs  compare  closely  with  those  of 
Jones  (1967)  for  propagation  below  an  isotropic  ionosphere  at  30  Hz. 

4.2  Cavity  Fields  from  Vertical  Electric  Sources 

Based  upon  the  discussion  in  Chapter  IH,  we  can  approximate  the 
lightning  stroke  as  a  vertical  electric  dipole  located  at  the  surface 
of  the  earth.  A  large  body  of  literature  on  radiation  from  dipoles  in 
the  earth-ionosphere  waveguide  has  been  published  and  only  the  salient 
features  of  the  approach  to  the  problem  will  be  presented  here.  For 
details  the  reader  is  referred  to  papers  by  Wait  (i960),  and  Galejs 
(1964a),  and  to  books  by  Wait  (1962)  and  Budden  (1961). 

Following  the  approach  of  Wait  (1962),  the  fields  can  be  obtained 
from  a  single  scalar  function,  If  the  cavity  has  conplete  -ynmetry 
in  the  0-direction  the  expressions  for  the  field  components  are 
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The  function  y  was  given  by  Wait  (1962) 

,  1M(»)  pv(-cose) 

4kha  sinvn 

(4.5) 
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The  symbols  in  expressions  (4.1)  -  (4.5)  have  the  following  meanings 
T,  =  characteristic  impedance  of  free  space  =  */\lq /cq 
k  =  free  space  propagation  constant  =  w  •fv-Qc0 
h  =  height  from  surface  of  earth  to  "lower  boundary"  of 
ionosphere 
a  *  radium  of  earth 
M(w)  s  source  dipole  moment 

6  =  angular  separation  between  source  and  receiver 
Pv(-cos6)  --  Legendere  function  of  complex  order,  v 

The  dimensionless  quantity  v,  which  is  a  function  of  frequency,  depends 
upon  the  properties  of  the  ionosphere  and  will  be  discussed  below.  In 
essence,  $  is  proportional  to  the  radial  component  of  the  Hertz  potential. 

Substituting  (4.5)  into  (4.1),  and  using  the  fact  that  Pv(-cos0)  must 
satisfy  the  Legendere  equation, 

§g  [sine  §e  pv(-cosG)J  +  v(v+i)pv(-cos0)  =  o 


we  obtain 
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(4.6) 


The  expressions  (4.5)  and  (4.6)  should  be  summed  over  all  allowed  values 
of  v,  corresponding  to  all  possible  waveguide  modes;  but  in  the  ELF  case, 


only  one  mode,  the  lowest  order  TM  mode  can  propagate  in  the  earth-iono¬ 
sphere  waveguide.  All  other  modes,  both  TM  and  TE,  are  well  below  their 
cutoff  frequencies  (Sudden,  19ol).  A  useful  representation,  called 
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Dougall's  expansion  (Magnus,  et  si.;  19&>)  of  the  Legendere  function  is 

m 

pv(-<!°8e)  _  i  V  a»i  .  ,  ... 

sinvB  R  L*  n(n+l)-v{  v+lj  Pn  00 

n=o 

where  the  Pn(cos0)  are  the  Legendere  polynomials,  so  that  the  radial  com¬ 
ponent  of  the  electric  field  becomes 
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(4.7) 


The  quantity  v  is  determined  by  the  conditions  on  E  and.  H  imposed 
by  Maxi-roll's  equations  at  the  earth  and  the  ionosphere.  Equivalent  to 
these  boundary  conditions  is  the  modal  resonance  equation  (Budden,  1961), 

R.(0)R  (0)e"2ikhcose  „  i  (4.8) 

— ■  & 


which  must  be  satisfied  if  a  waveguide  mode  is  to  be  self-consistent. 

In  (4.8),  0  is  the  complex  angle  of  incidence  (measured  from  the  normal 
to  the  reflecting  boundaries)  and  R_.  (0)  and  R  (9)  are  the  Fresnel  re- 

A  o 

flection  coefficients  at  the  ionosphere  and  earth,  respectively.  Wait 
(1962)  has  shown  that  v  +  =  ka  sin  0  so  that  if  allowed  values  of  cos  0 

are  determined  from  (4.8),  allowed  values  of  v  are  then  known.  It  should 
be  clear  that  v,  as  determined  by  (4.8)  is  a  function  of  frequency. 

Hence,  at  certain  frequencies  the  denominator  of  a  term  in  the  summation 
in  (4,7)  can  be  very  small  resulting  in  a  large  value  of  E  .  These  en¬ 
hancements  of  the  fields  at  certain  frequencies  were  first  predicted  by 
Schumann  (1952)  and  have  in  recent  years  been  called  "Schumann  resonances". 


-.1 


The  horizontal  electric  field  component,  E,,  is  very  small  since  at 
EL?  the  earth  is  a  very  good  conductor  and  the  electric  field  must  be 
very  nearly  ncmal  to  the  surface  of  the  earth.  For  this  reason,  the 
TM  mode  under  discussion  can  also  be  approximated  as  a  T01  mode  (Madden  - 
and  Thompson,  1965). 

The  "horizontal"  component  of  the  magnetic  field,  is  found  by 
expanding  the  potential  function,  f,  by  Dougall ’  s  representation  and 
substituting  the  result  into  (4.3); 


H  .1X1*1. 

0  4han 


2m- 1 

n(n+l)  -  v(v+l) 


iL_  p 

de  *n 


(cos0) 


mO 


(4.9) 


It  is  apparent  that  also  displays  the  resonant  behavior. 

The  preceding  discussion  can  be  summarized  as  follows;  El?  fields 
from  a  vertical  dipole  source  consist  essentially  of  a  vertical  E  field 
and  a  horizontal  H  field  At  the  earth  surface ,  the  magnetic  field 
lines  are  concentric  circles  about  the  source.  This  result  is,  however, 
only  true  if  the  cavity  formed  by  the  earth  and  the  ionosphere  is  per¬ 
fectly  symmetric  and  isotropic. 


4.3  Fields  From  A  Distribution  of  Independent  Sources 

The  results  of  the  preceding  section  will  now  be  generalized  to 
include  the  case  of  a  distribution  of  vertical  dipoles  on  the  surface 
of  the  earth  which  represents  a  reasonable  model  for  a  thunderstorm 
region.  Ac  this  print  we  will  no  longer  consider  the  field  components 
themselves,  but  th~  squares  of  the  magnitudes  of  the  field  components 
which  are  proportional  to  power.  The  relations  (4.7)  and  (4,9)  can  be 
thought  of  as  describing  linear  systems  where  M(oj)  is  the  input  signal 
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and  E  and  are  the  output  signals.  The  transfer  function  of  this  sya- 
r  v 

tem  depends  on  9,  the  angular  distance  between  source  and  receiver.  The 

case  of  several  lightning  strokes  can  then  be  thought  of  as  a  multiple 

input  linear  system  and  it  can  be  shown  (Appendix  B)  that  if  the  inputs 

to  such  a  system  are  uncorrelated,  the  total  power  at  the  output  is  the 

sum  of  the  powers  due  to  each  input.  Thus  if  it  is  ass, sued  that  the 

lightning  strokes  are  uncorrelated,  their  powers  can  be  added. 

The  square  of  the  absolute  value  of  the  E  field  due  to  a  distribu-  . 

2 

tion  of  dipoles  over  an  infinitesimal  area  a  sin0d0d0  located  at  a  dis¬ 
tance  9  from  the  receiver  on  the  earth  surfs'*.;  is  then  equal  to 


jEj2  -  «(«)  -vQr^ 


h 
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ttwO 


2w-l 

n(nfl)  -  v(v+l) 


Fn(coa0) 


a  sin0d0d0  (4.10) 


where  g(w)  is  the  dipole  moment  squared  per  unit  area,  which  is  in  the 
most  general  case  a  function  of  9  and  0.  In  a  more  abbreviated  form. 


|Ej2  -  GE(0,p)a2sin0d9# 


(4.11) 


For  further  development  we  assume  that  the  source  is  centered  at  an 
angular  position  0,  0  as  shown  on  figure  (4.1).  The  edge  of  the  source 
region  nearest  to  the  receiver  is  at  distance  (0-A)  and  the  far  edge 
at  (0+A).  The  azimuthal  spread  of  the  source  is  from  0-A^  to  0+^j. 

The  square  of  the  electric  field  at  the  receiver  due  to  this  source  is 
then  given  by 


«  0+A.  ©+A 

2  j.  0  .  P 

IeI  =  1  G„(0’,0)a  sin0'd0»d0'  (4.12) 

Va/0-a 


The  expression  for  the  magnetic  field  is  similar  but  the  H  field 


due  to  each  elauent  of  the  source  Is,  in  general,  in  a  different  direc 
tion-  If  we  let 


OHO,0) 


«(») 


2n+l 


Wx  L  n(ntl)  -  v(v+l)  n 


(cos0) 


n-cO 


<4.r 


then  the  square  of  the  magnetic  field  at  the  receiver  is  given  by 

Ih^I2  a  f  f  Ou(e',0*)  sin^ (0 * -Cf)a2sin9’ dS" d0'  (4.1J 

0  J0-^  Je-A  H 

2 

where  the  factor  sin  (0'-o)  accounts  for  the  vector  simulation.  The 
angle  0'  is  defined  in  Figure  (4.1)  so  that  the  source  region  extends 
in  the  azimuthal  direction  from  0'  *  0-A^  to  0*  »  0+&^  j  the  angle  or 
indicates  the  Inclination  of  the  receiving  coil  axis  with  respect  to 
the  0'  (and  0)  reference  line  which,  for  example,  could  be  the  east- 


west  direction. 


V.  TEE  RATIO  METHOD  FOR  SOURCE  LOCATION 


As  is  evident  from  the  discussion  of  the  preceding  chapter,  the 
received  Eli?  fields  are  dependent  upon  the  parameter  v,  which  is  de¬ 
termined  by  the  ionosphere,  and  the  parameter  0,  the  source-receiver 
separation.  In  the  source  location  problem,  it  is  advantageous  to  on- 
ploy  those  properties  of  the  received  signal  which  are  mainly  dependent 
on  8  and  are  relatively  independent  of  v  so  that  precise  knowledge  of 
the  state  of  the  ionosphere  is  not  required.  It  was  suggested  by 
Folk  (1969)  that  the  relative  sizes  of  adjacent  peaks  in  the  Schumann 
spectra  would  be  a  useful  Indicator  of  thunderstorm  location. 


5.1  The  Dominant  Ten  Approximation 

As  introduction  and  to  provide  insight  into  the  method,  we  Initially 

th 

make  a  radical  approximation.  Let  us  assume  that  at  the  n  resonant 
frequency  all  terms  except  the  nth  in  the  simulations  in  (4.7)  and  (4.9) 
are  negligible.  This  is  equivalent  to  assumi"*  that  the  earth-ionosphere 
cavity  is  very  nearly  lossless.  Under  this  approximation,  the  electric 
field  at  the  receiver  is 
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M(«n) 

4xa^w  e  h 
n  0 


v(v+l)(2n»l) 
n(n+l)  -  v(v+l) 


P  (eos9) 
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(5.1) 


and  the  magnetic  field  component  exciting  the  receiving  coil  is 


H_. 


lM(un) 

4xha 


(2m-l)sin(d-cr)  d  _ 
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at  the  n  resonant  frequency.  Prom  (5.1)  and  (5.2)  it  readily  follows 
that  at  the  n  resonant  frequency,  the  vertical  electric  field  varies 
as  P  (cos  0)  and  the  horizontal  magnetic  field  varies  as  P  (cos  0), 
This  behavior  is  illustrated  in  Fig.  (5*1)  for  the  first  two  resonances 
and  it  should  be  clear  from  Fig.  (5.1)  that  the  relative  sizes  of  the 
resonances  are  indicative  of  the  source-receiver  distance.  For  example, 
if  the  receiver  is  90°  from  the  source,  the  ratio  of  the  E  field  signal 
at  the  second  resonant  frequency  to  the  E  field  signal  at  the  first 
resonant  frequency  should  be  much  larger  than  the  corresponding  ratio 
if  the  receiver  is  65°  from  the  source.  For  the  readers'  convenience 
expressions  for  the  first  three  Legendere  polynomials  are  listed* 

P  ^  (cos  0)  -  cos  0 

P2  (cos  0)  -  \  (5  cosVl)  (5.3) 

(cos  0)  -  £  (5  cos50-3  cos  0) 

Equations  (5.1)  and  (5.2)  can  be  expressed  in  a  more  simple  form 
in  the  following  manner*  The  resonances  are  determined  by  the  condition 

v(v+l)  a:  n(n+l)  (5.^) 

and  since 

v  ♦  i  -  kaS  (5.5) 

the  resonant  frequencies  are  given  approximately  by 


(5.6) 


if  ReS  »  ImS.  If  the  cavity  has  perfectly  conducting  walls,  S  *  1, 
and  the  corresponding  resonant  frequencies  are 


(5.7) 
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so  that 


u>_ 


no 

ReS 


(5.8) 


Using  (5.4),  (5.5)  and  (5.8),  the  magnitude  of  the  eleotrlo  field  at  the 
th 

n  resonance  can  be  expressed  as 


M(.  )(2»1) 
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*A  Pn  ») 


(5.9) 
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The  frequency  dependence  of  S  for  several  different  ionosphere 
models  has  been  determined  by  Jones  (1967).  In  Figs.  (5.2)  and  (5.3) , 


a  -  -  0.182  f  ImS  d3/ton 


(5.10) 


and 


f-ReS  (5.11) 

where  f  is  the  frequency  in  Hertz,  c  is  the  velocity  of  light  in  air, 
and  v  is  the  phase  velocity  of  the  field  components  in  the  earth-iono¬ 
sphere  waveguide.  Note  that  if  the  cavity  walls  are  perfectly  conducting, 
the  attenuation,  a,  should  be  zero  which  is  confirmed  by  (5.10)  since 
ImS  «  0  for  the  perfect  cavity.  Fig.  (5.2)  reveals  an  interesting  prop¬ 
erty  of  the  earth-ionosphere  waveguide*  it  is  a  "slow  wave"  system  in 
the  Schunann  band,  i.e.  the  phase  velocity  is  always  less  than  c. 

The  first  three  3chumann  resonances  occur  at  about  8,  14,  and  20 
Hertz  with  approximately  6  Hertz  separating  them.  From  Figs.  (5.2)  and 
(5.3)  i*  is  evident  that  ReS  and  ImS  are  fairly  constant  in  the  frequency 
interval  between  adjacent  resonances.  The  ratio  of  the  electric  field 
magnitude  at  the  N  +  1  st  resonance  to  the  electric  field  magnitude  at 
the  nth  resonance  is  then  independent  of  the  properties  of  the  iono¬ 


sphere  and  is  given  by 


37 


Ifati  6) 

I  En  '  0) 


(5.12) 


Similarly,  it  is  not  difficult  to  show  that,  under  the  approximations 
that  have  been  made,  lE^/fcJ  is  also  independent  of  ionospheric  prop* 
erties  and  is  given  by 


I  3a  I  9>  (5W) 

i  Hn  1  'M(Vl)(^l)^li0^Pn(=o.  3) 

Another  parameter  which  will  be  of  interest  in  this  report  is  the 

XU  XU 

ratio  of  the  S  field  at  the  n  resonance  to  the  H  field  at  the  n 
resonance,  can  ue  ?«ai'iy  shown  to  be 
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(5.14) 


which  displays  a  dependence  on  the  ionosphere  model  (since  it  involves 
3)  but  has  the  interesting  feature  that  it  is  independent  of  the  light¬ 
ning  stroke  frequency  spectrum. 

The  preceding  discussion  applies  only  to  an  infinitesimal  point 

source  and  it  is  desirable  to  extend  these  results  to  a  distributed  source 

which  is  a  better  description  of  a  thunderstorm  region. 

Following  the  approach  in  Chapter  IV,  and  choosing  the  source  as 

in  Fig.  (4.1)  the  received  electric  and  magnetic  fields  in  a  high  Q 

cavity  can  be  determined  as  follows.  If  we  use  the  notation 
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UM)  -  *(«J  — |£§ —  P  (cos  e)  (5.16) 

no 

where,  as  before,  g(w'  Is  the  lightning  dipole  moment  squared  per  unit 
area;  the  magnitude  squared  of  the  received  electric  field  1& 

t+L0  0+fc 

iEn|2  -  J  J  M®’.  *’)a2  sin  0*d0*d0*  (5.1?) 

0-<d0  0-A 

and  the  magnitude  squared  of  the  magnetic  field  component  exciting  the 
receiving  coll  Is 

0+A0  0fA 

IhJ2  -  f  \  Gn(GS*')sin2(0*-or)^»0*d0'd0*  (5.18) 

n  0-A0  J0-A  H 


Assuming  that  g(u)  la  Independent  of  6  and  0,  the  squared  electric 
field  ratios  at  adjacent  resonances  In  the  high  Q  cavity  are 

v  2  w*  ,,22  P*  (eos0')sin0»d0* 

n  ^(w^Caw-l)  j1  Ap2(cos0»  )sin0»d0» 

0-A  n 

(5.19) 


and  the  squared  magnetic  field  ratios  are 


2  Q&  W<“*8’)J2  ^•d0’ 


^(Ujj)  (2n+l)  wrw.1>0  j  *jjL  Pn(cos0»)]  sinO’dO’ 
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The  impedeneea  o?  electric  field  to  magnetic  field  ratios  for  a  dis¬ 
tributed  source  in  the  high  Q  cavity  are  given  by 
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It  is  evident  from  equations  (5*19)  end  (5*20)  that  the  electric 

field  power  ratios  and  the  magnetic  field  power  ratios  are  independent 

of  the  azimuthal  spread,  A*,  of  the  source.  They  do  depend  on  the  dis- 

v 

tance,  9,  to  the  center  of  the  source  and  also  depend  on  the  range 
spread,  A,  however.  The  E  to  H  power  ratios,  eq.  (5.21),  depend  on 
A#,  A,  #  and  6  hut  are  not  a  function  of  the  lightning  spectrum  since 
they  are  ratios  of  quantities  measured  at  the  same  frequency.  In  the 
next  section,  we  dispense  with  the  assumption  of  a  high  Q,  cavity  and 
obtain  more  realistic  expressions  for  the  ratios. 

5.2  Ratios  for  Distributed  Sources  in  the  Lossy  Cavity 

Since  the  quality  factor  of  the  earth-ionc *^'iere  resonant  cavity 
is  of  the  order  of  4  (Jones,  1967)  the  assunpllon  that  only  the  reson¬ 
ance  tern  in  equations  (4.7)  or  (4.9)  is  significant  is  really  not  Jus¬ 
tified.  Relaxing  this  condition,  the  electric  field  ratios  are  given 
by  the  expression 
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and  the  ratios  of  the  magnetic  field  squared  at  the  n  +  1  st  resonance 
to  the  magnetic  field  squared  at  the  xrth  resonance  are 
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where  0^  (0,0)  is  given  by  eq.  (4,13).  Finally,  the  electric  to  mag¬ 


netic  field  ratios  are  given  by 
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Assuming  that  the  surface  distribution  of  dipoles  is  constant  over 
the  source;  or  equivalently,  that  g(u>)  does  not  depend  on  0  or  0,  it 
follows  that 
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There  are  four  important  parameters  to  be  determined  In  the  thunder¬ 
storm  location  problem  i  0,  the  angular  distance  or  range  to  the  center 
of  the  sourcef  0,  the  bearing  or  azimuth  of  the  source  center;  A»  the 
range  spread  of  source;  and  A0,  the  azimuthal  spread  of  the  source. 

The  E  ratios  and  H  ratios  contain  Information  about  6  and  A  but  are  in¬ 
dependent  of  0  and  A0.  The  E  to  H  ratios  depend  on  6,  A,  0,  and  A0. 

In  the  next  chapter  methods  are  discussed  for  estimating  these  parameters 
from  Schumann  resonance  data. 

Numerical  results  for  equations  (5.26)  -  (5.28)  can  be  obtained  by 
specifying  an  ionosphere  model  and  a  lightning  spectrum.  The  lightning 
spectrum  of  equation  (3.5)  has  been  used  successfully  by  Galejs  (1961a) 
and  will  be  employed  here. 

Jones  (1967)  has  shown  that  the  effect  of  heavy  ions  in  the  lower 
regions  (below  56  ton)  of  the  ionosphere  must  be  included  in  any  iono¬ 
sphere  model  to  be  used  below  50  Hertz,  Curves  ”0"  on  Pigs.  (5.2)  and 
(5.3)  give  ReS  and  ImS  for  an  ionosphere  profile  determined  from  reaction 
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rates  of  atmospheric  constitrmnt  gases  by  fierce  a sd  Cole  (1965).  This 
profile  includes  the  effect  of  the  heavy  ions  and  values  of  v  can  be 
computed  from  the  curves  "e*  by  v  a  kaS  -  0.5* 

The  variation  with  source  distance  of  the  ratios  of  the  adjacent 
resonant  peaks  in  the  received  electric  field  power  spectrum  is  shown 
in  Pig.  (5*4)  for  A  »  20°.  The  results  of  calculations  using  both  the 
hlgh-Q,  dominant  tens  approximation  and  the  full  series,  inhomogeneous 
ionosphere  model  are  shown.  Note  that  the  relative  sizes  of  the  peaks 
display  a  marked  variation  with  source-receiver  distance.  It  is  also 
evident  that  the  ratios  for  the  high  ft  cavity  are  not  radically  dif¬ 
ferent  from  the  ratios  for  the  realistic  cavity.  It  is  also  apparent 
from  Pig.  5. 1»  that  the  results  for  the  point  source  and  distributed 
source  are  qualitatively  the  same,  l.e.  at  8  -  90°  the  second  resonance 
is  large  compared  to  the  first  resonance,  etc.  For  this  reason  and 
because  the  choice  of  the  ionosphere  model  is  not  extremely  critical, 
it  is  evident  that  eqs.  (5.12,  5.1?»  5.14)  can  yield  great  insight  into 
the  source  location  problem. 

The  magnetic  field  power  spectrum  ratios  for  both  the  high  ft  and 
low  ft  cases  are  shown  in  Fig.  (5.5).  Here  the  ionosphere  model  seems 
to  have  a  greater  effect  but  qualitatively  the  results  are  again  the 
same  for  both  ionospheres. 

The  expression  for  the  electric  to  magnetic  ratios,  eq.  (5.28) 
is  fairly  complicated  in  that  it  depends  on  all  four  of  the  source 
parameters.  For  the  moment  we  simplify  the  expression  by  choosing  the 
special  case  A0  small,  a  -  0°,  and  0  ■  90°.  This  effectively  assumes 


that  the  azimuthal  spread  is  small  (leas  than  15°)  and  that  the  receiving 
coil  is  oriented  so  that  its  maximum  response  is  in  the  direction  of  the 
center  of  the  source.  Under  these  assumptions  the  coefficient  multiply¬ 
ing  the  ratio  of  integrals  in  (5.28)  is  unity.  For  this  special  case 
the  ratios  \?  /Hnj  are  plotted  as  a  function  of  0  for  both  the  high  Q 
and  low  Q  cavity  in  Fig.  (5.6).  Again  the  qualitative  results  are  the 
same  for  both  ionosphere  models. 

5.3  Effect  of  Ionosphere  Model  on  Ratios 

It  is  evident  from  the  results  of  the  previous  section  that  even 
a  very  drastic  change  in  the  properties  of  the  ionosphere  such  as  the 
difference  between  the  Pierce-Cole  profile  and  the  perfectly  conducting 
model  does  not  fundamentally  alter  the  variation  of  the  ratios  with 
source-receiver  distance.  To  investigate  this  effect  in  a  more  meaning¬ 
ful  manner,  ratio  calculations  baaed  upon  two  different  inhomogeneous 
ionosphere  profiles  were  performed. 

The  exponential  conductivity  profile  used  by  Galejs  (l96la)  and 
the  Pierce-Cole  (1965)  conductivity  profile— both  for  day  ionospheres— 
are  shown  in  Fig.  (5.7).  The  Galejs  profile  is  based  upon  particle 
density  and  collision  frequency  data  for  electrons  only  (Galejs,  196lb)j 
the  Pierce  and  Cole  curve,  on  the  other  hand,  includes  the  effects  of 
positive  and  negative  ions  which  account  for  the  fact  that  the  two 
models  have  conductivities  that  differ  by  as  much  as  two  orders  of 
magnitude  at  the  lower  altitudes, 

A  comparison  of  the  ratios  for  the  Pierce-Cole  (or  Jones)  iono¬ 
sphere,  the  Galejs  ionosphere  and  the  perfectly  conducting  ionosphere 
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is  given  by  Pigs.  (5.8  -  5>.l*).  A a  would  be  expected,  the  ratios  for 
the  Galejs  model  are  closer  to  the  results  of  the  Piero  e-Cole  (Jones) 
model  than  to  those  for  the  high  Q  oavity.  The  perfectly  conducting 
case  isf  of  course,  a  very  unrealistic  approximation  to  the  ionosphere 
and  should  only  be  used  to  provide  initial  insight  into  the  Schumann 
phenomena. 

If  an  experimental  spectrum  is  given  and  the  ratios  are  computed 
from  it,  the  source  distances  estimated  from  the  Jones  or  Gala J a  models 
would  differ  at  most  by  about  10°  which  is  an  acceptable  error  for  a 
source  of  A  -  20°  (total  source  width  4C°). 

The  relative  Invariance  of  the  ratios  of  spectral  peaks  with  iono¬ 
sphere  properties  indicates  that  observations  of  these  ratios  might  be 
practically  useful  for  determining  source-receiver  separation.  It  has 
been  pointed  out  by  Large  and  Wait  (1968)  and  by  Balser  and  Wagner  (1962) 
that  the  values  of  the  Schumann  resonant  frequencies  themselves  depend 
also  upon  source-receiver  separation}  it  should  be  noted  however,  that 
the  exact  values  of  the  resonant  frequencies  are  very  strongly  dependent 
on  the  properties  of  the  ionosphere.  Illustrative  of  this  fact  are 
Pigs.  (5.15  -  5»20)  where  the  six  resonant  frequencies  are  shown  as  a 
function  of  source  distance  for  the  ionosphere  models  of  Galejs  and 
Pierce-Cole. 

Suppose  one  were  to  attempt  to  determine  source  location  from 
resonant  frequency  measurements:  if  the  first  resonance  in  the  magnetic 
field  were  at  8,0  Hz  (for  example),  the  estimated  source  distance  would 
be  6j5°  based  upon  the  Pierce  and  Cole  model  and  110°  based  upon  the 
Galejs  model,  a  difference  of  47°  (5000  km).  If  the  second  resonance 
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in  the  H  field  were  at  13.7  Hz,  the  distances  would  be  28°  and  75° »  also 
a  large  variation.  Note  that  for  a  higher  or  lower  peak  frequency  only- 
one  of  the  models  could  yield  a  source  distance;  in  this  case,  we  could 
use  the  resonant  frequency  to  choose  the  ,fbest"  of  the  two  ionosphere 
profiles.  The  ionospheric  dependence  of  the  peak  frequencies  of  the 
second  and  third  resonances  is  strikingly  illustrated  by  Pigs.  (5.17), 
(5.19),  ard  (5.20)  which  could  also  be  used  to  determine  which  of  the 
two  ionosphere  profiles  is  "best". 

It  is  clear  that  the  resonant  frequencies  are  determined  at  least 
equally  by  the  ionosphere  profile  and  the  source  location.  In  fact, 
for  the  higher  resonances,  the  ionosphere  seems  to  be  the  dominant  factor 
in  the  resonant  frequency  variations.  To  obtain  sourc e-receiver  informa¬ 
tion  from  resonant  frequency  measurements  would  require  detailed  know¬ 
ledge  of  the  local  as  well  as  global  properties  of  the  ionosphere;  the 
ratio  method,  however,  demands  only  fairly  crude  ionospheric  data. 


VI.  PROCEDURES  FOR  THE  LOCATION  OP  SOURCES  BY  THE  RftTIO  METHOD 

Having  demonstrated  that  the  ratios  of  the  spectral  peak  powers 

are  relatively  independent  of  the  properties  of  the  ionosphere,  we  now 

turn  to  a  discussion  of  methods  for  estimating  the  source  parameters 

3.  0.  A.  and  A.  from  measured  ratio  data. 

0 

According  to  Jones  (1969)  a  given  spectrum  corresponds  to  one 
and  only  one  source-receiver  separatior,  i.e.  the  spectra  are  unique. 

It  is  reasonable,  therefore,  to  assume  v.hat  a  given  ordered  sequence 
of  measured  ratios  (|Eg/1SjJ%  |EyE2|2,  (Hg/HjJ2,  |Hy'H2|‘)  can  corres¬ 
pond  to  only  one  source  distance,  6.  Consequently,  one  would  expect 
that  unique  determinations  of  source  receiver  separation  from  ELF 
spectra  is  possible. 

6.1.  General  Discussion  of  Source  Distance  Estimation 

Suppose  we  wish  to  find  the  "beat  match  distance"  between  n  ratios; 
R^»  R^  ...  Rn.  For  example,  R1  might  be  \'R2/Ej\2,  Rg  »  |Hyi^|2,  R^  « 

| Eg/Ei |  ,  etc.  It  is  desired  to  find  one  angle,  8,  which  is  the  solution 
to  the  n  equations, 

\  =  fx(e) 

R2  -  f2(e) 

r  =  f  (e) 

n  nv  ' 


(6.1) 


4? 


The  f  functions  are  basically  the  'ratios  of  squared  suns  of  Legendere 

o 

polynomials  for  or  r^*08  °f  squared  sums  of  deriva- 

2 

tives  of  Legendere  polynomials  for  j  hence,  the  functions 

^(0)  are  non-linear  and  single  valued,  but  the  inverse  functions 
^“"(R^)  are  non-linsar  and  multivalued. 

In  the  discussions  below,  let  0^  be  the  i^  root  for  the  ratio 
and  let  Nj  be  the  nuriber  of  roots  of  Rj  *  fj(0),  *,e*  3  at 

0  *  G^,  0i>  ...  0^  .  If  there  are  n  ratios  to  be  matched,  the  solu- 
tiona  of  equations  (6.1)  will  yield  then  sets  of  roots 

0, 1 ,  Gg1,  O^1,  ...  G^1  by  solving  R1  =  ^(0) 

0^,  022,  0?2,  ...  0^2  b/  solving  Rg  =  f?(0) 


e^.  02n»  •••  ®n  n  fcy  solving  R^  =  fn(e) 

n 

a  total  of  (N^  +  N,,  +  ...  Nn)  angles  of  which  only  one  from  each  set 
is  the  desired  an«?.e.  Our  task  is  to  choose  the  n  values  of  0^  that 
are  mutually  closest  remembering  that  ..nere  is  some  measurement  error 
in  the  given  values  of  By  mutually  closest  is  meant  those  n  values 
of  0^  for  which  pj  and  qn  are  selected  such  that 


«PJt  Qjj)  = 


n  n 


j=l  n>k>j 


(6.2) 


is  a  minimum  1  £  Pj  £  Nj 


1<  ^<Nq 
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The  distance  estimate,  6,  is  than  specified  by 

e  »-  (el  +  e2  ♦  ...  e11)  (6.3) 

n 

where  the  selected  0^  are  the  n  mutually  closest  roots. 

The  ratio  values  will,  of  course,  not  be  exact  due  to  the  noise 
like  properties  of  ELP  signals.  Is  is  discussed  in  Chapter  U,  the 
power  spectrum  estimates  are  random  variables  with  standard  deviation 

where  p  is  the  expected  value  of  the  power  spectrum  estimates.  For  a 
frequency  resolution  of  0.3  Be  and  a  time  sample  13  nilnutes  long  the 
standard  deviation  is  about  five  percent  of  the  expected  value.  The 
probability  that  a  normally  distributed  random  variable  is  less  than  2 
standard  deviations  from  its  expected  value  is  about  .95.  Therefore, 
it  is  quite  certain  that  the  actual  ratios  are  within  plus  or  minus  ten 
percent  of  the  estimated  ratios. 

If  the  possibility  of  errors  in  ratio  estimates  is  included  in  the 

source  location  problem  we  define  a  perturbed  set  of  ratios  given  by 

(l  +  (l  +  ^^^2'  •••  (l  +  En)Rn  where  the  are  in  the  range 

-.15  <  +.15  .  The  method  is  then  to  determine  Q,  for  all  different 

1.  2  3  xi 

combinations  of  these  perturbations.  The  roots  0  ,  9  ,  0  ,  ...  0 
which  give  the  smallest  ft  are  then  used  to  determine  the  best  match 
distance  by  (6.3). 

The  larger  perturbations  are  less  probable  than  the  smaller  pertur¬ 
bations,  i.e,  it  is  most  likely  that  the  measured  ratios  are  equal  to 
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the  actual  ratio.  Because  of  this,  It  Is  desirable  to  employ  a  "stop¬ 
ping  rule"  as  follows.  Rather  than  searching  for  the  absolute  wrtmimai 

Q  as  above,  we  stop  searching  for  a  match  when  the  maximum  of  \B'l  -  0*  | 

J  oqJ 

for  all  j  and  k  is  less  than  some  small  specified  angular  distance  (5 
for  the  A  *  20°  model,  for  exanple).  The  search  for  the  best  match  is 
done  beginning  with  the  smallest  perturbations  on  the  roots  and  working 
towards  the  larger  perturbations.  The  ordering  of  the  perturbation  se¬ 
quence  is  determined  by  the  quantity 


Fki 


(6.4) 


i-1 


so  that  the  perturbation  vectors  (elt  e^,  ...  eq)  which  give  the  smallest 
p  are  used  first  in  the  root  search  procedure,  then  the  perturbation  set 
which  gives  the  next  larger  y  is  used  and  so  on,  until  either  a  match 
is  achieved  or  all  possible  perturbations  (within  the  specified  limits) 
have  been  investigated.  This  algorithm  must,  of  course,  be  performed 
on  a  digital  computer. 


6.2.  The  Distance  Estimation  Procedure  for  Matching  Two  Ratios 

As  a  simple  application  of  the  general  method  outlined  above  we  con¬ 
sider  the  problem  of  finding  the  unique  6  which  gives  the  best  match  to 
two  given  ratios  and  R^.  R^  could  be  )  E^/E^  | 2  and  Rg  could  be  |  Eyl52 1  * 
in  the  A  =  20°  model,  for  example.  The  procedure  requires  frequent 
evaluation  of  the  inverse  functions  f^”1^);  i  »  1,  2,  ...  n  (n  ■  2  In 
this  example)  which  is  done  by  linear  interpolation  on  tables  of  values 
of  0  and  R^  and  tables  of  values  of  0  and  Rg.  The  multivalued  property 
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of  the  functions  ^(0)  and  fg(0)  can  be  handled  by  dividing  tte  tables 
into  several  piecewise  singlevalued  regions  of  6.  In  the  ample  under 
discussion  f^O)  -  jS^/S^I2  and  f2(0)  -  |E/S2|2  where  fj(0)  has  four 
mono  tonic  regions*  0  <  55°»  55°  <  0  <  90°,  90°  <  0  <  125°*  and  0  >  125°* 
The  function  f2(0)  also  has  the  same  four  mono  tonic  regions* 

0  <  55°,  55°  <  0  <  90°,  90°  <  0  <  125°,  and  0  >  125°.  The  reader  can 
verify  these  statements  by  referlng  to  Figure  (5.4).  Hence,  for  a 
given  and  Rg,  there  will  be  at  most  four  roots  of  ^(0)*  0^,  02, 

9y  and  four  roots  of  fg(0)*  6^*  ©2t  ®4* 

The  root  search  procedure  is  as  follows.  The  smallest  perturbation 
pair  elt  e2  is  chosen  (these  will  both  be  zero)  and  the  roots  in  each 
monotonic  region  of  each  ratio  are  found  by  interpolation  on  the  tables 

of  theoretical  ratio  calculations.  Then,  the  smallest  of  the  quantities 

12  o 

1©^  -  0j  |  is  found  and  this  difference  compared  to  5°.  If  it  is  less 

than  5°  it  is  termed  a  match  and  the  distance  estimate,  0  is  determined 

by 

1 

0  -  \  (@"  +  02) 

1  2 

where  0  and  0  are  the  two  closest  angles,  the  first  selected  from  the 

1111  2  2  2 
set  (©!  t  ©2  »  »  ®4  )  and  the  second  from  the  set  (0^  ,  02  ,  0,  ', 

0^  ).  Of  course,  the  closest  two  angles  may  be  more  than  5  apart. 

If  this  is  the  case  the  next  larger  perturbation  on  and  Rg  is  used 
to  determine  the  roots,  the  difference  between  the  closest  two  roots 
is  compared  to  5°  and  so  on,  until  either  there  is  a  match  or  all  per¬ 
turbation  pairs  have  been  tried.  In  the  latter  case,  the  model  is 
changed  (for  example  a  different  &  is  used)  and  the  procedure  is 
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initiated  again.  1  flowchart  of  the  method  for  Batching  two  ratios  is 
shown  in  Pigure  (6.1).  Refering  to  the  flowchart,  MAXKPS  is  the  largest 
perturbation  considered,  COPS  is  the  inc resent  between  different  per¬ 
turbations,  KEHO  is  the  maxlnnm  allowed  distance  difference  for  a 
match,  KPT  is  the  magnitude  of  the  angular  distance  between  the  two 
closest  roots  for  a  given  permutation  pair,  and  NPERTS  is  the  total 
nunfeer  of  different  perturbation  pairs. 

6.3.  Estimation  of  Source  Distance  and  A  from  Ratio  Data 

The  source  range  spread  can  be  estimated  from  the  E  ratios  and/or 
H  ratios  in  the  following  way.  Suppose  it  is  desired  to  choose  the  best 
distance,  6,  and  best  match  A  from  the  possible  models  A  *•  0°,  6  -  5°» 

A  «*  10°,  A  «  20°.  The  procedure  is  merely  an  extension  cf  the  method 
described  in  section  2i  one  determines  the  best  0  from  the  4«0° 
model  and  stores  the  perturbations  required  for  a  match.  Then  the  same 
is  done  for  the  other  three  A  values.  Whichever  A  value  gives  a  match 
with  the  smallest  perturbation  is  the  best  match  A  and  the  corresponding 
9  is  the  best  match  distance. 

6.4.  Determination  of  Azimuth  and  Azimuthal  Spread  of  the  Source 
Refering  to  figure  (4.1)  let  the  0*0  line  from  the  receiving 

site  be  the  east-west  line.  A  coil  whose  axis  is  along  the  east-west 
line  (cr  «  0)  will  be  excited  by  the  magnetic  field  component  given  by 

|H^|2  -  [A0  -  £  cos20sin2A0]  3Cn(9,A)  (6.3) 

9fA 

Kn(0,A)  Gfi  (e»)a2sin0«de* 


where 


and  the  magnitude  squared  of  the  magnetic  field  component  measured  by 
a  coil  whose  axis  is  along  the  nortb-south  line  (or  «  9C  )  is  given  by 

JH*13!  -  [A#  ♦  £  cos20sin2A0]  Xn(0,A)  (6.6) 

The  electric  field  magnitude  squared  at  the  receiver  is,  as  before. 


Dividing  equation  (6*5)  by  equation  (6.6)  one  obtains 


!  H® 

n 

Mm 

A 

0  -  &  cos20sin2A^ 

H 

NS 

n 

*  A 

.  +  A  cos20sin2A. 
#  *  0 

and  if  is  small  (A^  <  10°)  this  reduces  to 


(6.8) 


(6.9) 


A  plot  of  against  0  for  small  is  shown  in  Pigure  (6.2). 

Note  that  for  a  fixed  |hPV  ^|2  value  there  are  four  possible  aai- 
asuthal  locations,  0.  The  number  of  possible  source  bearings  can  be 
reduced  to  two  as  follows.  Suppose  two  coils  with  perpendicular  axes 
are  mounted  on  a  turntable  so  the  angle  or  between  the  coils  and  the 
0  »  0  line,  see  Pigure  (4.1),  can  be  varied.  Let  jH^J  be  the  magni¬ 
tude  squared  of  the  magnetic  field  component  exciting  a  coil  whose 

A 

axis  makes  an  angle  a  ■  or^  with  0*0  line  and  let  |  be  the  magnitude 
squared  of  the  magnetic  field  component  exciting  a  coil  whose  axis  is  at 
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an  angle  cr2  ■  cr^  +  */2  with  respect  to  the  0  «  0  line.  Then  it  is 
evident  from  the  discussion  in  Chapter  5  that 


A.-  \  cos2(0-ar)sin2A  2 

A*+  1"  cos2(0-«)sin2A^  ~  ^ 


(6.10) 


if  A .  is  small.  Hence,  at  or  -  0  (or  0  +  x)  the  ratio  of  the  fields 
0 

measured  by  the  two  colls  must  be  small  and  by  rotating  the  colls  until 

2 

if  is  a  minimum,  the  axis  of  coil  1  lies  on  the  great  circle  path 

between  source  and  receiver. 

Hie  assumption  of  small  A .  is  not  critical.  Even  in  the  highly 

0 

unlikely  case  b.  »  30°.  0  estimates  based  on  the  A.  *  0°  model  will  be 
0  0 

in  error  by  only  about  5°  or  10°  as  is  shown  in  Figure  (6.3)  where 

Eq.  (6.10)  for  various  A*  is  plotted. 

0 

We  have  seen  that  estimates  of  9  and  A  can  be  obtained  from  the 
ratios  | Kg/fcJ2,  \Kf&2\2  and/or  (H^/Sj2,  iHy^j2  aid  that  by  employ¬ 
ing  a  pair  of  orthogonal  coils  the  azimuth  oj.  the  source  can  be'  found. 

2  2 

From  either  of  the  quantities  |En/H^|  or  IVfl  the  value  of  A^ 
can  be  obtained.  For  example. 


E 


n 

a® 

n 


2A 


2  ■  _  n 


c  (e,A) 


A^-4  cos2^8in2Aii  TC  (0,A) 
0  c  0  n 


(6.11) 


where  curves  of  £n(0, A)/Kh(0,A)  versus  0  for  various  A  are  given  in 
Chapter  5.  After  having  found  0,  A,  and  0  by  the  methods  described 


above  eq.  (6.11)  contains  only  one  unknown,  A^,  which  can  then  be 


determined  if  measurements  of  are  available  (although  it  is 

clear  that  (6.11)  shows  only  a  very  small  dependence  on  Atf). 
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Thus,  all  the  significant  parameters  of  a  thunderstorm  source  can 
be  determined  by  measurements  at  a  single  receiving  station.  There  is, 
however,  an  ambiguity  of  l80°  in  the  determination  of  the  source  azimuth 
which  could  probably  be  resolved  by  the  elimination  of  unlikely  (e.g. 
temperate  zone  winter)  source  locations. 


6.5.  An  Alternate  Method  for  Source  Distance  Estimation  from  Measured 
Ratios 

The  solution  of  the  problem  poced  by  equations  (6.1)  can  also  be 
treated  as  a  minimization  of  the  distance  between  two  vectors  in  an  n 
dimensional  vector  space.  Let  the  vector  of  given,  experimental  ratio 
values  be  denoted  by 


(6.12) 


and  the  vector  of  function  values,  which  depend  on  both  0  and  A  as 
parameters,  be  denoted  by 


fx(efA) 

f2(e,A) 


fni(M) 


(6.15) 


The  distance  between  these  two  vectors  is  given  by  the  expression 


■sswwHW1 
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.-/It*  -*Mrf 


and  this  is  a  ndnimun  for  the  same  9  and  A  that  minimize 


Q  *  2<  K  -  fi(0^)32 


(6.14) 


The  solution  of  equations  (6.1)  is  then  where  Q.  »  0.  There  are,  of 

course,  some  errors  in  the  measurement  of  the  quantities  r., ,  r ^  ...  r 

l  c.  n 

so  that  there  may  not  exist  any  (6,4)  pair  for  which  Q  =  0.  One  would 
expect,  however,  that  useful  estimates  of  the  source  distance  and 
range  extent,  A,  could  be  obtained  by  minimizing  Q. 

The  function  Q  in  equation  (6.l4)  has  the  undesirable  property 
that  a  small  error  in  a  large  ratio  value  can  be  an  equal  or  even 
larger  contributor  to  the  total  mean  squared  error  than  a  large  error 
in  a  small  ratio  estimate.  It  would  be  preferable  to  weight  all  the 
errors  on  a  percentage  basis  which  can  be  accomplished  bTr  employing  a 
minimizing  function  of  the  form 


S  rr,  -  f,(e,A)-,2 

*  -  L  Lr^r - J 


(6.13) 


The  standard  method  for  determining  the  9  and  A  which  minimize  the  ex¬ 
pression  in  equation  (6.15)  would  be  to  solve  the  two  simultaneous  non¬ 
linear  equations 


22-0 
se  '  0 


If-o 
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i’or  8  and  A.  The  expressions  for  the  f,(0,A)  are,  however,  extrenely 
involved  and  the  standard  approach  seems  unwarranted.  The  method  used 
here  is  to  simply  guess  the  8  and  A  which  minimize  Q.  The  procedure 
would  be  as  follows:  choose  a  value  of  A,  say  0°.  'Then  evaluate  Q 
for  every  value  of  8  from  0°  to  l8o°  in  steps  of  one  degree.  The  value 
of  0  for  which  Q  is  a  minimum  would  be  the  best  0  for  that  particular  A. 
Next  choose  another  value  of  A,  say  5°,  and  repeat  the  above  search 
for  the  minimum  Q.  The  procedure  is  continued  until  all  values  of  A 
have  been  investigated  and  their  corresponding  best  match  0's  have  been 
found.  The  solution  would  be  that  (8, A)  pair  for  which  Q  is  a  minimum. 
The  algorithm  outlined  above  will  always  find  a  solution  to  R  =  P(0) 
regardless  of  how  poor  a  match  is  obtained.  A  rule  to  determine  whether 
or  not  a  solution  has  been  found  would  be  to  test  the  quantity  \/Q^~/n 
(average  normalized  deviation)  ’gainst  a  prescribed  tolerance,  say  .05. 

6.6.  Estimation  of  Source  Distance  Prom  |E  /H  j  ratios 

1  n'  n1 

The  ratio  of  the  electric  field  power  spectral  density  to  the 
magnetic  field  power  spectral  density  at  their  respective  nth  resonances 
is  essentially  independent  of  the  lightning  spectrum  since  the  resonances 
of  the  two  fields  occur  at  very  nearly  the  same  frequency.  A  measurable 
quantity  (which  is  not  the  total  magnetic  field  magnitude  squared)  is 
determined  by  the  addition  of  eqs.  (6.5)  and  (6.6)  to  be 
EW  ?  NS 

|H  i  +  |H  r  =  M  (e,A) 

:i  n 
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and  the  electric  field  magnitude  squared  Is 


so  that 


\\\  -  ^0  *„<«•*> 


i-fi2  ♦  hf  i2 


^(e,A) 


2  2  2  2 

Hence,  the  ratios  [&J  2/(l^"i  ♦  |H®|  )i  |E2|2/(|^M|  ♦  |H®|  )  and 

|E//(|Hf|2  *  |«f  I2  ')  depend  only  on  the  source  distance  and  distance 

spread,  not  on  the  lightning  spectrum,  the  source  azimuth  or  azimuthal 

spread.  These  remaining  source  parameters  can  be  determined  as  described 

in  section  (6.4)  and  it  is  clear  that  all  the  parameters  of  the  source 

can  be  estimated  regardless  of  the  lightning  spectrum  provided  that 

ratios  taken  at  the  same  frequency  such  as  |En/Hn|  ,  1^  I  * 

lEX“l2-  ”  ivf  I2  are  used.  This  suggests  a  method  of  estimating 

the  lightning  spectrum  at  the  Schumann  resonant  frequencies  since  the 
2  2 

ratios  |E  and  do  depend  on  the  lightning  frequency 


n+r  n* 


spectrum. 


VII.  LOCATION  OP  TWO  ELF  SOURCE  REGIONS  PROM  A 


59 


Since  there  are  now  eight  parameters  which  describe  the  vocation  and 
extent  of  the  two  sources,  it  should  be  expected  that  estimation  of  these 
parameters  would  be  more  difficult  than  in  the  previous  single  source 
case. 


7.2  Estimation  of  Range-Parameters  for  the  Two  Source  Model 

It  is  advantageous  to  employ  in  the  location  procedure  those  quan¬ 
tities  which  depend  on  the  least  number  of  independent  parameters  since 
their  determination  will  be  then  greatly  simplified.  The  electric  field 
spectrmn  depends  upon  6  of  the  8  source  parameters  while  each  of  the  two 
magnetic  field  spectra  is  a  function  of  the  entire  set  of  8  parameters . 
This  complexity  can  be  reduced  by  combining  equations  (1.2)  and  (7-3)  t° 
obtain 

|n“|2  ♦  !h“3|2  -  24, *  (ef  .  4a)  +  24^  M  (8b,  4g)  (7.4) 


The  number  of  variables  can  be  further  reduced  if  ratios  of  + 

and/or  |e|2  are  considered.  For  example, 


2\  t,.  *  \ 


2d .  ^  +  2k 

0A  n+l  PB  n 


e* 


(7.5) 


where  =  e  (©A*  k^)  (see  eq.  6.6)  is  evaluated  at  the  (rw-l)st  reson¬ 
ance  of  the  total  received  electric  field  (since  the  total  electric  field 
is  what  can  be  measured).  Equation  (7.5)  can  be  written  as 


pA 

'"'n+l 


6  e; 


B 


n+l 


£*  + 
n 


6 


B 


n 


(7.6) 


rO 


%*here  6  =  / 1 ).  The  ratios  in  (7.5)  thus  depend  on  9^,  0^,  Ag, 

B  A 

and  6, 

We  can  further  reduce  the  number  of  variables  to  three  by  making  the 
assumption,  AA  =  Ag  =  A  *  5°  corresponding  to  a  great  circle  distance  2A» 
of  about  700  miles,  which  is  probably  a  representative  or  average  value 
for  large  thunderstorm  areas  (Kealohs,  et.  al.,  1958;  Byers,  19^9).  In 
any  event,  the  ratios  display  only  a  second  order  (compared  to  0) 
variation  with  A  in  the  range  0°  <_  A  <,  10°  which  probably  includes  all 
realistic  source  sizes.  As  can  be  seen  from  Figures  (7*1)  and  (7.2), 
a  A  value  of  5°  is  representative  for  the  range  0°  <,  A  <.  10°. 

The  following  seven  ratios  are  then  dependent  only  on  0^,  0^  and  6» 


(7.7) 


(7.8) 


n  n 

/  +  6  H0  ' 

n  n 


n  ■  1, 


2,  3 


(7.9) 


The  above  seven  ratios  are  not  all  independent;  the  particular  subset 
to  be  used  will  depend  on  the  available  data.  Also,  the  ratios  of  the 
set  (7.9)  are  independent  of  the  lightning  spectrum  (provided  the 
source  spectrum  at  location  A  is  the  same  as  at  location  B),  If  the 
selected  ratios  are  denoted  by  f. (©.,  0_,  6),  i  =  1,  2,  ...n,  the 
parameters  9^,  0g,  and  6  could  be  estimated  by  minimizing  the  quantity 


>#•»  *  w r**>  *v^>‘ 
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n  2 

ft  *  7  Oi  (C±  ♦  6  D1)  -  (^  +  6  B±)J  (7.14) 

i=l 

which  leads  to  the  following  linear  equation  for  6 
n  n 

T  (C±  *±  -  A±)  (Di  r±  -  B1)  4  6  7  (D±  r±  -  B±)2  «  0  (7.15) 

i=l  i*l 

so  that 

n 

l  (A,  -  C,  ri)(D1  r,  - 

6  =.  iji -  (7.16) 

Z  <Di  ri  -  vs 

i=l 

since 

S  -  Z  <di  ri  -  v2  *  °  ^ 

d8  i-i 

we  are  guaranteed  that  this  6  minimizes  ft  for  a  specified  0A  and  6^. 

The  method  of  determining  0^,  8g,  and  6  from  measured  ratios  can  now 
be  outlined  as  follows: 

We  investigate  all  values  of  8^  and  6g  between  10°  and  100°  in 
steps  of  a  prescribed  increment,  say  2°.  The  first  value  of  0^,  say 
10°,  is  chosen  and  then  a  value  of  0g.  All  the  quantities  on  the  right 
hand  side  of  (?.l6)  are  then  specified  and  a  value  of  6  can  be  calculated. 
The  corresponding  value  of  ft  is  determined  from  (7.14)  and  stored.  Then 
tha  next  @n  is  chosen,  6  determined,  and  ft  calculated.  If  this  ft  is 

IJ 

less  than  the  previous  ft,  the  former  (along  with  the  associated  values 
of  0  ,  0g,  and  6)  is  stored  in  place  of  the  latter  and  the  next  Qg  in 
sequence  is  used  to  determine  a  ft.  This  procedure  is  continued  until 


«*»•**&>*■ 


all  possible  values  of  6g  are  investigated.  The  next  value  of  6^ 
is  chosen  and  all  possible  values  of  8g  are  used  to  determine  a  Q 
value  and  so  on  until  the  entire  allowed  set  of  9^,  0g  pairs  have 
been  used  to  evaluate  the  6  which  yields  a  minimum  Q.  The  values  of 
0^,  9g,  and  6  which  produce  the  ndnimun  of  this  set  of  minimums  are  the 
best  choices  of  these  source  parameters  for  the  given  ratios. 

7.3  Estimation  of  the  Azimuthal  Parameters  in  the  Two  -  Source  Model 
There  are  fifteen  ratios  of  spectral  peak  powers  which  depend  upon 
the  azimuths  or  bearings  of  the  two  sources;  they  are: 


|  Mf 


;  n  *  1,  2,  3 


|v?1  1  ”  •  3 

Mff  *  n  =  1,  2,  3 

KlAltff * Ml2) * 


!  n  *  1,  2,  3 


MW  *  Ml2) 


I  n  *  1,  2,  3 


where  n  denotes  the  ordinal  number  of  the  resonance.  After  examining 
these  quantities,  however  it  becomes  clear  that  a  knowledge  of  some 
of  them  implies  a  knowledge  of  the  others.  In  fact,  values  for  all 
of  them  are  specified  by  measurements  of  any  two  members  of  the  set 


|«f.  W  KM 


we  shall  confine  our  attention  to  the  first  two  ratios  of  this  set; 
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a  different  choice  will  not  affect  the  augments  which  follow. 

By  applying  equations  (7.1  -  7*3)  these  ratios  are 

2  ^ASin  ^  ]  >£  *  \  J  "d 

L\4=o»  aVA  J  **  ♦  [y4««  a^sin  24^  ]  )£ 


» 

L 

5®! 

n  1 

se 

L 

(7.13) 


E 


H®| 
n  1 


2a.  e4  +  2a.  e® 

0A  n  0g  n 


[\^°S  ^A^  \  j  +  [\^C°3  ^BSin  %  ]  *n 


(7.19) 


n  =  1,  2,  3  where 


*A,f* 

E4’8  *  Gu  (G»)sin0’d6* 

n  .  .  A  n 

0A.B~a  n 


(7.20) 


e4’8  =  f 

n  J 


\b*l 

g_  (e'jaii.e’de' 

eA.B-*  “ 


(7-21) 


and  expressions  for  G^  and  Gg  were  given  previously  (equations  (4.13)  and 
(5.23). 

A  word  about  the  geometry  used  in  the  thunderstorm  location  pro¬ 
cedure  is  in  order.  Since  the  area  of  a  thunderstorm  region  which  is 

2 

2A  in  azimuthal  spread  and  2A  in  range  spread  is  4a  A  sin  9  sin  A 
v  $ 

where  a  is  the  radius  of  the  earth  (6.4  x  10'  km),  the  area  of  the  source 


depends  upon  the  source-receiver  separation  9.  Hence,  at  very  close 
(near  0°)  or  very  distant  (near  l80°)  stoims  the  source  area  for  a  fixed 
A  becomes  unrealistically  small.  This  can  be  accounted  for,  however, 


** 


f 

•  I 
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by  allowing  A.  to  assume  fairly  large  values,  say  as  much  as  40°.  These 
v 

remarks  become  Important  in  determining  the  validity  of  certain  approxi¬ 
mations  which  are  to  be  discussed  below. 

The  equations  (7.18)  and  (7.19)  can  be  greatly  simplified  by  making 

the  approximation  that  A^  and  A^  are  small  enough  so  that  the  factors 

A  B 

sin  2  Am  and  sin  2  AM  can  be  replaced  by  2A  and  2A*  ;  this  applies 
9A  eA 

if  A.  and  A^  are  not  larger  tnan  20°.  We  can  preserve  a  representative 
®A  ®B 

source  area  of  the  order  of  at  most  300,000  square  miles  (roughly  the 
size  of  the  state  of  Texas)  in  the  range  10°  £  0  £  YJ0°  by  adjusting  the 
A^*s  within  the  upper  limit  of  20°. 

The  simplified  forms  of  eqs.  (7.18)  and  (7.19)  are 


H® 

n 


sin\  +  6  **  sln20B 


n 


cos20a  +  6  eos20 


(7.22) 


E 

n 

n 


+ 


6  6® 
"n  n 


JjJ  sin20A  +  6  «? 


sin20n 
n  H3 


(7.23) 


where 


n  =  1,  2,  3 


6  *  K  A. 

®B  »A 


Since  we  have  previously  shown  how  to  estimate  ©A»  0^  and  6  the  only 
unknowns  in  (7.22)  and  (7.23)  are  0A  and  0g,  provided  measurements  of  the 
quantities  or  the  left  hand  side  are  available.  Subsequent  discussion 
will  require  some  algebraic  manipulation  and  we  shall  simplify  the  notation; 
let  = 


r2  = 


h/V3 


r3  " 


H f/Ef 


Kl  = 


E, 


A 


EW 


FW  **  EW  ^ 

K2  =  jVV  .  X3  »  Ey'Rj  . 


Using  the  identity  sin  x  +  cos  x  =  1  the  equations  (7.22)  and  (7*25)  are 
seen  to  be  six  simultaneous  equations  which  are  linear  in  the  two  un¬ 
knowns  sin20,  and  sin^Z*  .  The  first  three  of  these  equations  corresponding 
to  (7.22)  are 


»n  *in\  *  5  *n  sln\  ■  7TI  (*£  *  5  $ 

n 


(7.24) 


n  =  1,2,5  and  the  three  equations  corresponding  to  (7*25)  are 


*n  ”ln\  «  5  »n  ,in\  =  K~(4  + 6  en5 

n 


n  =  1,2,5. 


(7.25) 


It  is  clear  that  the  equations  (7*24)  and  (7*25)  are  identical  since 
it  follows  from 


n  /  A  ,  Bn  1  /_A  ,  .Bv  in 

r  +1  ^n  +  6  *n^  =  K  ^n  +  6  ~  2Aa 

n  n 


(7.26) 


hence,  we  need  deal  with  only  one  of  the  two  sets  of  equations,  either 
(7.24)  or  (7.25),  to  determine  0^,  and  0g.  Without  loss  of  generality, 
we  shall  choose  to  solve  equations  (7.24)  and  shall  express  them  in 
the  matrix  form 


MS  =  C 


(7.27) 


where 


M  « 


(7.28) 
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and 


8  = 


sin^d. 


sin 


L 

r. 


(7.29) 


C  * 


rl+l 


(Hi4  +  6  »!B) 


2  /*.  A  p  y  B\ 

r^+l  ^2  +  ®  ^2  ' 


r,+l  (*3A  +  6  M3B) 


(7.30) 


As  is  well  known  (Henmerle ,  1967),  the  best  solution  to  (7.27)  in  the 
least  squares  'tense  is  obtained  by  solving  the  equation 


jfus.lfc  (7.31) 

where  the  superscript  T  denotes  the  transpose  of  the  matrix.  The 

matrix  equation  (7*31)  represents  two  simultaneous  equations  in  the 
2  2 

two  unknowns  sin  0  and  sin  0g  and  can  be  solved  by  conventional  means. 
It  is  inst inactive ,  however,  to  examine  the  determinant  of  the  matrix 
of  coefficients,  mF  M,  in  (7. 31).  This  determinant  is  given  by  the 
expression 


det  (wF  M)  =  62  {  [(a-jV  +  (*2A)2  +  (*3A)2  K)2  +  (»2B)2  +  (*3B)2  j 


-  (X\B  +  *2\B  +  »„V  ]  } 


3 

(7.32) 


A  B 

which  vanishes  if  0  =  9  .  Therefore  if  the  two  source  centers  are  at 


equal  distances  from  the  receiver,  ratios  of  the  type  shown  in  (7.22)  or 
(7.23)  can  provide  no  information  on  the  azimuthal  location  of  the 
sources.  This  fact  becomes  clear  if  one  examines  the  expressions  for 


the  ratios  given  by  (7*22)  or  (7*23)  under  the  condition  that  0A  *  0g 
(or  *4  "  *n  “  *n  anct  Jn  *  Cn  *  ffn^  *ince  the  r*tio8  become  independent 
of  the  resonance  number.  In  factf  if  ©A  «  00  and  AA  at  Ag 


I*/  .  240  e*  *  a,  fn  -  »*,, 

A  ** 

2 

if  -  “».5lAA^  +  ^  »^#n  *  24 '«  *iA'»n 
1  A  B 

2 

-  2^  cos20aK^  +  2A*  cos^^j®  «  cos^'jin 
A  B 


where 


V  *  ^  “  ^0 

A  B 

and  0'  ie  an  equivalent  -  fictitious  -  azimuth  for  the  two  sources  located 

at  0.  =  0n  with  0.  0_. 

A  B  A  B 

A  single  received  spectrum  is  produced  which  is  indistinguishable  from  the 
spectrum  from  one  source. 

7.^  Detecting  the  Presence  of  More  Than  One  Major  Source 

If  a  single  major  source  is  exciting  the  cavity  the  ratios 

^  -£  coa2<6e in2&0 
J^S  "  ^  +£  cos20sin2A0 


are  obviously  independent  of  the  n,  the  ordinal  number  of  the  resonance. 


For  the  case  of  two  sources,  however, 


|«|  given  by  equation 
(7*18)  are  clearly  dependent  on  n  unless  =  6g  and  2=  Ag.  We  can 
then  make  the  following  statement  1  if  the  measured  quantities 
are  different  for  different  resonance  numbers,  there  must  be  more  than 
one  source.  It  does  not  follow,  however,  that  if  these  ratios  are 
independent  of  n  there  is  only  one  source,  since  as  we  demonstrated  in 
the  previous  section,  two  sources  at  equal  0  but  different  azimuth  0 


can  imitate  one  source  at  a  single  receiving  station 


VUI.  EFFECTS  OP  BACKGROUND  NOISE  ON  THE  SOURCE 
LOCATION  METHOD 

Up  to  this  point,  it  has  been  assumed  that  the  measured  ELF  spectra 
are  not  contaminated  by  signals  of  an  origin  other  than  major  intense 
thunderstorm  regions.  This  is,  of  course,  not  always  true  and  in  this 
chapter  the  effects  of  various  types  of  noise  on  the  estimated  source 
distances  will  be  discussed. 

8.1  Effects  of  Extra-Terrestrial  Noise  on  Source  Location  Estimates. 

Schumann  spectra  c-f  the  magnetic  field  ere  frequently  observed  to 
have  a  sharply  increasing  power  level  with  decreasing  frequency.  This 
effect  can  be  pronounced  enough  tc  cause  errors  in  the  estimates  of  the 
power  levels  at  the  lower  Schumann  resonant  frequencies  (and  consequent 
errors  in  the  estimation  of  source  location)  and  can  even  be  so  strong 
as  to  entirely  remove  any  evidence  of  resonances  (see  Figure  8.1).  The 
source  of  this  noise  is  probably  outside  the  earth-ionosphere  cavity  and 
possibly  involves  hydromagnetic  emission  phenomena;  it  cannot  be  pro¬ 
duced  by  distant  thunderstorms  since  the  analytical  expression  for  the 
magnetic  field  does  not  contain  a  "zeroth"  resonance  term.  In  contrast, 
the  electric  field  does  have  a  "zeroth"  resonance. 

Numerical  experiments  were  conducted  to  determine  the  effect  cf 
noise  of  this  type  on  the  location  estimates.  The  frequency  spectrum 
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of  the  noise  can  be  approximated  by 

H(f)  =  Af"®  (3.1) 

where  or  is  in  the  range  1  to  2.  It  is  clear  that  for  a  given  A,  the 
larger  a  is  the  less  a  Schumann  spectrum  will  be  contaminated.  If  the 
spectrua  resulting  from  thunderstorm  activity  is  designated  by  T(f), 
the  total  received  spectrum,  S(f),  (including  the  noise)  can  be  written  as 

S(f)  *  N(f)  +  T(f)  («:?} 


provided  N  and  T  are  the  power  spectra  of  uncorrelated  random  processes. 
Assuming  a  source  location  and  strength,  one  can  then  evaluate  S(f)  for 
different  values  of  A  and  or  and  use  these  “noisy"  spectra  in  the  source 
location  procedure.  Figures  (3.2)  and  (8.3)  illustrate  the  results  of 
those  investigations.  In  these  graphs  8q  is  the  true  location  of  the 
thunderstorm  source;  the  abcissa  is  the  ratio  cf  the  spectral  power  at 
2  Hz  to  the  spectral  power  at  the  first  resonant  frequency  i.e. 


where  f^  is  the  apparent  first  resonant  frequency.  Since  thunderstorm 
produced  power  at  2  Hz  is  negligible  compared  to  the  noise  power  at  2  Hz, 


(8.4) 


C  was  chosen  in  this  form  so  that  it  could  be  used  as  a  criterion  for 
accepting  or  rejecting  a  spectrum  as  useful  for  the  location  of 
thunderstorms , 

Variation  along  the  abscissa  can  also  be  roughly  interpreted  as 
variation  of  the  percent  noise  power  at  the  apparent  first  resonance. 
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The  reason  for  this  is  as  follows:  The  percent  ncise  power,  P,  at  the 
first  resonance  ia  given  by 

H(0 

p  -  100  £  K<fx> ;  £  (8-5) 

and  since  is  approximately  a  constant  (f^  depends  on  A,  a,  and  the 
source  location  and  strength;  but  the  range  of  variation  is  at  most 
about  +_  1  Hz) ,  we  have 

P  =  100  [jf(f1)A(2)]C  (8.6) 

so  that  P  and  C  are  related  by  a  constant  if  a  is  fixed.  We  cannot  ex¬ 
perimentally  observe  P,  but  we  can  measure  C.  Figure  (8.2)  shows  the 
variation  of  the  estimated  source  distances  9  for  a  =  1  and  varying  C 
or  P  (which  corresponds  to  varying  A).  The  estimated  source  locations 
are  plotted  until  the  estimates  are  more  than  10°  in  error.  The  amount 
of  noise  which  is  tolerable  is  dependent  on  0Q.  For  example,  if  the 
source  is  at  120°  distance,  source  distance  estimates  are  unreliable 
if  C  >  .1  (P  >  12  V°)/but  if  the  source  is  at  20°,  C  <  4  or  P  nearly 
100  °/o  does  not  affect  the  distance  estimates. 

Selection  of  Of  =  1  is  probably  a  too  severe  assumption  about  the 

noise  frequency  variation.  Experimental  data  indicate  that  more  often 

2 

a  a  2.  Figure  (8.5)  describes  the  effect  of  noise  if  N(f)  -  A/f  ; 

0^  =  120’  is  still  the  worst  case  but  not  until  C  _>  4  (P  >,  25  °/°)  must 
we  reject  a  spectrum. 

It  can  be  concluded  that  the  estimation  procedure  is  fairly  insen¬ 
sitive  to  this  type  of  noise. 
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8.2  Uniform  Equatorial  Distribution  of  Low  Level  Thunderstorm  Activity. 

In  this  section,  we  examine  the  effect  on  the  source  location  pro¬ 
cedure  of  low  level  thunderstorm  activity  distributed  over  a  large  equa¬ 
torial  area.  According  to  Pierce  (1969a),  large  scale  thunderstorm  areas 

have  a  characteristic  dimension  of  about  1000  Kra  (A  =  4.5°)  or  an  area  of 
6  2 

10°  tan.  Herman  (1963)  has  concluded  that  the  world’s  major  thunderstorm 

6  2 

regions  have  an  area  of  5  x  10°  tan"  (about  l/500th  of  the  earth’s  surface 
area).  These  thunderstorm  regions  consist  of  complexes  30-40  tan  in  dia¬ 
meter  and  separated  by  about  100  km  (lierce,  1969b).  On  the  average,  a 
specific  complex  will  contain  one  or  two  active  cells  which  are  6  tan  in 
diameter  and  produce  three  lightning  strokes  per  minute  (Pierce,  1969b). 
Using  Pierce’s  model  we  can  estimate  the  flash  rate  of  a  typical 

major  thunderstorm  area.  Imagine  that  a  typical  complex  is  placed  at 

4  2 

the  corners  of  a  square  region  100  tan  x  100  tan.  This  10  tan  region 
then  contains  one  complex  and  on  the  average  will  produce  6  flashes  per 
minute  so  that  one  of  our  typical  major  thunderstorm  regions  will  pro¬ 
duce  600  flashes/minute.  Assuming  that  two  of  these  major  .regions  are 
active  at  a  given  time,  the  flash  rate  from  these  two  major  regions  would 
then  be  1200  discharges  per  minute  or  20  discharges  per  second  from  the 
world  major  thunderstorm  regions.  It  is  generally  believed  that  there 
are  three  major  thunderstorm  regions— tropical  America.  Africa,  and 
Asia;  however,  meteorological  conditions  which  are  favorable  to  thunder¬ 
storms  depend  on  local  time  (late  afternoon  is  usually  considered  opti¬ 
mum  for  the  occurrence  of  thunderstorms)  so  that  at  most  two  major  re¬ 
gions  are  lively  to  be  active  at  a  given  time  (Smith,  lQ6l) . 
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According  to  the  widely  accepted  thunderstorm  theory  of  atmospheric 
electricity,  woi ..d-wide  lightning  activity  is  the  generator  of  the  atmos¬ 
pheric  potential  gradient  (Chalmers,  196?)  and  from  measurements  of  this 
potential  gradient,  the  total  world-wide  lightning  activity  at  a  given 
time  can  be  estimated.  In  this  way,  it  can  be  roughly  determined  that 
there  are  at  least  100  .lightning  discharges/second  occurring  over  the 
entire  world  (Smith,  1961).  Polk  (1969) ,  using  a  related  technique,  lias 
estimated  that  there  must  be  250  '’typical'’  discnarges  per  second  over 
the  earth. 

Prom  these  very  rough  estimates  it  seems  that  the  amount  of 
lightning  activity  concentrated  in  the  major  thunderstorm  regions  is  only 
about  twenty  percent  of  the  total  world-wide  activity.  Suppose  it  is 
assumed  that  the  other  80  percent  (about  80  discharges  per  second)  is 
uniformly  distributed  in  an  equatorially  centered  band  between  50°  north 
latitude  and  ;0°  south  latitude.  Then  the  flash  rate  surface  density 
in  this  band  (its  area  is  exactly  one-half  the  total  surface  area  of  the 
earth;  about  2.58  x  10®  km?)  is  3.1  x  10-^  f lashes/km^/sec .  Comparing 

this  with  the  flash  rate  density  in  the  major  regions,  2  x  10*"^  flashes/ 

2 

km  /sec.,  we  see  that  the  source  strength  (g(w))  in  equation  (4.10)  in 
the  major  regions  is  about  two  orders  of  magnitude  larger  than  that  in 
the  uniform  equatorial  band.  It  remains  to  compare  the  received  power 
levels  from  the  uniform  band  of  scattered  "background"  sources  and  the 
major  regions. 

Suppose  there  are  strokes  uniformly  distributed  over  an  area  A. 
Then  the  total  dipole  moment  squared  over  this  area  is  N^M^(w)  and  we 


(3.7) 


N^fa)  =  j  ^  g(«)dA 


p  r 

1 

*  k 

and  If  g(w)  is  uniform  over  A, 

N.M2  («) 

*<«)  =  -T - 


(8.8) 


but  «  pA  where  p  is  the  stroke  surface  density  (strokes  per  unit 
area)  so  that 

sM  *  pM^{«)  (8.9) 

We  can  new  modify  equation  (3-3)  to  account  for  the  stroke  surface 
density  and  write 


g(»)  *  (const.)  0  exp(-9,l  x  lCf'  u»)  (8.10) 

where  the  constant  depends  upon  the  properties  of  our  "typical"  light¬ 
ning  stroke.  Knowledge  of  an  absolute  value  for  this  constant  is  not 
necessary  in  the  ratio  method;  only  the  frequency  variation  is  impor¬ 
tant  and, as  was  discussed  in  Chapter  HI,  the  frequency  variation  is 
held  in  fair  agreement  by  several  researchers. 

It  is  instructive,  however,  to  examine  the  likely  range  of  values 
of  g(u>)  so  that  calculated  signal  levels  can  be  compared  with  the 
measured  levels.  As  mentioned  in  Chapter  III,  the  average  return  stroke 
lias  a  peak  current  of  about  p0,000  amps  and  a  duration  cf  about  100 
micro-sec  nds  (Pierce,  1963;  Watt,  I960);  assuming  a  stroke  length  of 
two  kilometers,  the  peak  value  of  the  current  moment  of  the  return  stroke 
is  6  x  10 ^  amp-meters.  The  bandwidth  of  the  pulse  i3  of  the  order  of 
the  reciprocal  -f  the  pulse  length  or  10  KHz.  Therefore,  in  the  region 
0  <  f  <  ;0  Rz.  the  return  stroke  spectrum  can  be  considered  nearly  white 
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and  of  magnitude  equal  to  (peak  amplitude)  x  (pulse  length)  or  6  x  1GV 
amp -meter  secs.  Consultation  of  any  of  the  standard  tables  on  Fourier 
transform  pairs  show  that  the  detailed  pulse  shape  is  unimportant 
(e.g.  see  Korn  and  Korn,  1968)  for  the  purposes  of  estimating  expected 
signal  levels.  Note  that  this  value  of  6  z  10^  for  M(w),  although  ob¬ 
tained  in  a  very  rough  way,  compares  quite  favorably  with  the  value, 

5  x  105,  given  by  Pierce  (1963),  We  then  can  conclude  that  M^(w)  for 

6  7  2  2  2 

the  return  stroke  must  be  in  the  range  10°  -  1C  amp  -meter  -sec  so 
that  g(«)  would  be  as  follows: 

g(o)  for  uniform  lew  level  background  is  [j.l  x  10  ,  J.l  x  10~f] 
g(u>)  for  major  thunderstorm  areas  is  [2  x  10~\  2  x  10”^]. 

Using  the  largest  g(w)  in  their  respective  ranges,  a  major  source 

(A  =  5°,  A^  5°)  at  a  distance  of  40°  should  produce  theoretical  first 

resonance  field  amplitudes  as 

jE^j  =  2Ag£^(^0,5)  =  l8p  volts/meter 

i&J  =  2A^N^(40,5)  =  .03u  aznps/meter  (8.11) 

(or  .O^  milligamraa) 

o 

and  the  uniform  +_  O  band  of  low  density  "background"  thunderstorm 
activity  produces  fields  (at  Kingston,  R.  I.;  located  at  approximately 
4l.6°N,  71.8°W)  of  amplitudes 

JE, i  =  19. volts/meter 

JrLj  =  .0^5u  &  ps/meter 
(or  .056  milligamraa) 


(8.12) 
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These  values  are  one  or  two  orders  of  magnitude  smaller  than  those  re¬ 
ported  in  the  literature  (Polk,  19 69).  A  possible  reason  for  this  dis¬ 
crepancy  is  that  the  assumed  value  of  K(isj)  is  one  or  two  orders  of  mag¬ 
nitude  toe  small.  Rycroft  (1965)  has  obtained  good  agreement  between 
theoretical  and  experimental ly  measured  8  Ez  field  amplitudes  vising  the 
dipole  moment  M(w)  =  1,4  x  IQ"*  at  8  Hz  which  is  roughly  two  orders  of 
magnitude  larger  than  would  be  expected  for  an  "average"  return  stroke. 
Prom  "slow  tail"  measurements,  Hughes  (1967)  has  estimated  M(w)  to  be 
about  3  x  1(P  at  100  Hz  and  notes  that  this  is  about  40  dB  larger  than 
the  estimates  based  upon  the  return  stroke  only  (Williams,  1959).  In 
a  more  recent  paper,  Hughes  (1969)  provides  experimental  measurements  of 
M(ui)  which  if  extrapolated  down  in  frequency  to  the  region  of  the  first 
three  resonances  yield  a  value  of  about  10^  amp-meters.  Polk  (1969)  has 
used  estimates  of  the  vertical  electrostatic  field  to  conclude  that 
M(ui)/h  =  4':  amps  where  h  is  the  ionospheric  height.  For  h  =  50  km,  this 
yields  a  moment  of  M(u>)  =  2.2  x  10^  amp-meters.  The  empirical  formula 
of  Pierce  (1963)  describing  time  variations  of  the  different  components 
of  a  typical  discharge  (leader,  return,  and  continuing  currents)  can  be 
used  to  estimate  the  frequency  dependence  of  the  lightning  stroke  dipole 
moment.  As  was  shown  in  Figure  (3.1).  the  total  M^(w)  is  about  two  or¬ 
ders  of  magnitude  higher  than  that  for  only  the  return  stroke.  These 
discussions  seem  to  indicate  that,  as  was  pointed  out  by  Pierce  (1963). 
all  the  processes  of  a  lightning  discharge  are  important  in  the  excita¬ 
tion  of  the  Schumann  resonances  and  not  exclusively  the  return  stroke  as 
is  widely  assumed.  If  the  previously  used  value  of  M(u>)  is  increased  by 


a  factor  of  10,  the  theoretical  signal  levels  (equations  (8*11)  and 
(8.12))  compare  well  with  published  data  (Polk,  1969). 

Equations  (8.11)  and  (8.12)  indicate  that  the  signal  levels  from  a 
major  thunderstorm  area  and  from  the  equatorial  region  of  background 
noise  could  be  of  comparable  size.  This  raises  the  question  as  to 
whether  the  "signal1*  from  a  major  region  can  be  sufficiently  discrimin¬ 
ated  from  the  "noise"  of  the  equatorial  background  ac  that  the  location 
of  the  major  region  can  be  determined. 

It  is  reasonable  to  expect  that  the  background  thunderstorm  noise 
is  seasonally  dependent,  i.e,  the  background  noise  spectrins  would  be  dif¬ 
ferent  for  spring-fall,  northern  simmer  and  northern  winter.  Estimates 
of  the  geographical  distribution  of  the  equatorial  noise  belt  can  be 
obtained  from  the  literature  (Smith,  1961).  Assuming  a  distribution 
uniform  in  longitude,  the  extent  in  latitude  of  the  noise  belt  for  the 
various  seasons  is  roughly 

spring-fall  50°S  to  50°N 

northern  sunmer  10°S  to  60°N 

northern  winter  40°3  to  10°N 

The  electric  and  magnetic  field  power  spectra  received  in  Rhode  Island 
for  these  three  cases  are  shown  In  Figures  (8.4)  and  (0.5).  The  maxi¬ 
mum  power  level  variation  for  -che  different  seasons  is  only  about  thirty 

percent;  this  might  be  expected  since  the  source  area  for  the  three  cases 
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is  fairly  constant:  winter  area  =  2.1  x  10  km  ,  spring  fall  area  = 

8  2  8  2 

2.58  x  10'"  km  .  summer  area  =  2.68  x  10"  km  .  A  most  important  charac¬ 
teristic  of  these  spectra  which  should  be  pointed  out  is  the  fact  that 


the  peak  powers  at  the  resonances  decrease  with  increasing  resonance 
number,  i.e.  the  ratios  jE^/E^J2,  lEy^j2,  J H^/H^ | 2 »  jHyftgj2  are 
all  less  than  unity. 

One  may  ask  if  the  "real  world*  Schumann  spectra  ever  depart  from 
the  estimates  given  by  Figures  (8.4)  and  (8.5)  or  if  the  spectra  show 
a  day  to  day  variation  during  the  same  season.  Obviously,  if  the 
spectra  are  the  same  from  day  to  day  or  from  week  to  week,  the  location 
of  our  hypothesized  "major"  regions  would  be  impossible  because  either 
these  major  regions  do  not  exist,  or  their  level  is  so  low  that  their 
presence  cannot  be  detected  in  the  background  noise.  Actually  measured 
spectra  do  display  significant  day  to  day  and  even  hour  to  hour  variations. 

Typical  power  spectra,  measured  in  Rhode  Island,  are  plotted  in 
Figures  (8.6)  -  (8.8).  Figure  (8.6)  shows  that  in  the  space  of  six 
hours  a  significant  change  in  the  east-west  magnetic  field  power  spec¬ 
trum  occurred  on  February  19,  1967.  Both  the  ratios  and  the  resonant 
frequencies  show  marked  changes.  Figure  (8.7)  shows  a  similar  variation. 
An  outstanding  example  of  spectral  deviation  from  the  background  esti¬ 
mates  of  Figure  (3.4)  is  shown  in  Figure  (8.8)  where  the  electric  field 
power  spectrum  at  1^00  EST  on  August  14,  1968  is  plotted.  Here  the 
power  level  increases  with  the  order  of  the  resonance  and  it  is  clear 
that  at  this  time,  the  background  level  must  have  been  low.  Spectra  of 
this  type  are  not  at  all  rare  occurrences  and  we  can  conclude  that 
usually  the  background  noise  does  not  completely  conceal  the  signal 
from  the  major  regions. 

Equations  (8.11)  and  (8.12)  indicate  that  the  "noise"  level  pro¬ 
duced  by  the  equatorial  background  and  that  from  a  major  source  are  com- 
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parable j  however,  experiment  shows  that  significant  departures  from  a 
constant  ambient  spectrum  dc  occur.  In  the  paragraphs  to  follow,  we  shall 
examine  the  problem  of  how  much  noise  is  tolerable  in  the  location  of  the 
principal  thunders  tom  regions. 

iT  p 

To  simplify  the  analysis,  the  case  of  locating  a  single  10°  km 
major  source  in  the  presence  of  the  uniform  equatorial  background  region 
will  be  considered.  The  stroke  densities  in  the  major  and  background 
regions  are  denoted  by  and  Pgj  the  number  of  strokes  per  unit  time 
in  the  regions  are  and  Ng.  We  shall  use  Smith  *3  estimate  of  =  100 

discharges  per  second  over  the  entire  world  and  write 


or 


100  11  PMAM  +  PB*B 


(8.13) 

(8.14) 


where  and  Ag  are  the  areas  of  the  major  and  background  regions.  For 
the  purposes  of  computation 


and 


efi(w)  =  PgM^(ui) 


(8.15) 

(8.16) 


If  the  relative  flash  density  Pq/°^*  and  the  areas  A^  and  Ag  (dependent 
ot  the  season)  are  specified  the  p„  and  pn  can  be  determined  from  (8.l4) 
and  the  spectra  from  the  major  source  and  the  background  can  be  calcu¬ 
lated.  Figures  (8.9),  (8.10),  (8.11)  show  the  effect  of  the  background 
noise  on  the  source  distance  estimates  for  northern  winter,  spring-fall, 
and  northern  summer.  On  these  curves  the  abscissa  is  the  relative 


strength  of  the  background  source  Pg/P-,,  the  rdinate  is  the  estimated 
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source  distance  and  the  curves  are  parametric  in  the  actual  source  loca¬ 
tion,  8q.  If  Ajj  and  Ag  are  specified,  ^  *  PBAB  can  be  determined  for  a 
given  value  of  Vb  and  the  values  of  Ng  (strokes  per  3econd  in  the  back¬ 
ground  region)  are  displayed  at  the  top  of  the  page  for  the  regions  indi¬ 
cated.  The  curves  simply  display  the  deviation  of  the  estimated  location 
from  the  true  location  with  increasing  noise.  It  can  be  seen  that  there 
is  very  little  seasonal  dependence  of  noise  effects,  a  result  which 

might  be  expected  since  the  corresponding  noise  spectra  are  very  similar. 

-4 

We  can  conclude  that  Pg/Pjj  must  be  less  than  or  equal  to  about  2  x  10 
_<  5)  for  the  determination  of  precise  source  locations.  This  result 
is  consistent  with  the  calculations  of  (8.11)  and  (8.12)  where  it  was 

found  that  the  "noise"  level  and  the  "signal"  level  were  about  equal  for 

—2 

the  case  Pg/PM  —  10  .  This  means  the  received  signals  are  about  fifty 

percent  noise  and  it  is  not  surprising  that  the  resulting  estimated 
source  locations  are  very  different  from  the  true  locations.  This  re¬ 
sult  should  be  contrasted  with  the  case  of  the  Af~°  type  additive  noise 
whore  almost  50  percent  noise  at  the  first  resonance  was  tolerable?  but 
this  is  due  to  the  fact  that  the  exponential  noise  contaminates  only  the 
first  and  second  resonances  of  the  magnetic  field.  Both  electric  and 
magnetic  field  ratios  are  used  in  the  estimation  procedure.  On  the 
other  hand  equatorial  background  noise  contaminates  all  the  peak  powers 
of  both  the  electric  and  magnetic  field.  It  appears  that  the  ratio 


method  is  quite  sensitive  to  our  hypothesized  equatorial  ncise  sources. 

Figures  (8.6)  to  (8.8)  indicate,  however,  that  the  estimates  used 
in  equations  (8.11)  and  (8.12)  for  the  relative  strengths  of  uniformly 
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distributed  lightning  activity  as  compai*ed  with  activity  in  the  major 
thunderstorm  regions  are  probably  erroneous .  In  particular,  these  esti¬ 
mates  are  very  sensitive  to  the  area  of  the  major  thunderstorm  region; 
for  example,  if  the  assigned  area  of  a  major  region  were  increased  to 
5  x  10  km  (Herman,  1968)  the  "signal"  to  noise  ratio  would  be  greatly 
improved.  The  physical  parameter's  describing  thunderstorms  and  lightning 
have  large  variability  and  are  still  the  subject  of  active  research.  Con¬ 
sequently,  the  numerical  values  quoted  here  are  at  best  crude  estimates 
and  are  only  indicative  of  the  order  of  magnitude  of  the  particular 
quantity.  It  is  clear  that  the  assumed  uniform  equatorial  background 
noise  must  te  either  more  restricted  in  area  or  must  have  a  lower  rela¬ 
tive  flash  density  than  assumed  in  the  discussion  leading  to  (8.11)  and 
(8.12)  if  Elf  spectra  lead  to  the  correct  determination  of  thunderstorm 
regions.  To  the  extent  to  which  it  becomes  possible  therefore  to  con¬ 
firm  Elf  results  by  direct  surface  or  satellite  observations,  information 
is  also  obtained  about  the  relative  magnitude  of  "major''  and  "background" 
lightning  activity.  For  example,  if  it  should  turn  out  that  the  Elf 
method  gives  consistently  correct  results  during  certain  periods,  but 
not  daring  others,  it  would  be  possible  to  draw  conclusions  conclusions 
concerning  the  distribution  in  time  of  the  "background"  activity. 


.  LOCATION  OF  THUNDERSTORMS  BY  ANALYSIS  OF  SCHUMANN’  SPECTRA 


In  this  chapter  the  ratio  method  is  applied  to  observed  spectra  and 
estimated  locations  are  given.  Where  possible,  the  resulting  location 
estimates  are  compared  with  available  meteorological  data. 


9.1  Location  of  a  Single  Active  Thunderstorm  Region 

During  the  period  5-12  June,  196?  an  exceptional  amount  of  thunder¬ 
storm  aeti\  j.ty  oc cured  near  Salt  Lake  City  (Lao  h0°  kp'N,  Long  131°  52’ W) 
in  the  Rocky  Mountain  region  of  the  United  States  and  detailed  records 
of  the  associated  phenomena  have  been  published  (Kealoha,  et  al.;  1966), 

The  power  spectra  of  the  vertical  electric  field  and  the  east-west  magnetic 
field  received  in  central  Rhode  Island  (71°  4l'  01"W,  4lc  31*  53"N)  at 
17';'0  EST  to  17^2  EST  on  June  6,  196?  are  shown  in  Figure  (9.1).  This 
particular  time  was  selected  for  analysis  because  simultaneous  satellite 


cloud  cover  pictures  and  radar  precipitation  echo  data  were  available 
(Kealoha,  et  al.;  1968).  The  radar  data  indicate  that  at  this  time 


the  center  of  the  precipitation  region  for  this  thunderstorm  area  was 

0  'O  * 

at  distance  S  =  29  ,  and  bearing  <t>  -  -  lb  s,9ee  Figure  4.1  for  definition 
of  angles)  and  had  an  area  of  about  .8  x  10’  krr."  roughly  corresponding 
to  A 0  -  p . 5 " ,  A  =  ph.  Employing  the  minimization  method  described  in 
section  (6.5)  and  using  the  ratios  computed  frem  Figure  (9*1) ,  source 


locations  (6)  for  several  different  source  widths  (A)  were  estimated. 
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The  source-  distance  was  estimated  fx'Cff  the  ratios  ,  \E^/Epl'  ; 

jHfrf/}^j%  and  The  source  azimuths  were  obtained  by  re¬ 

arranging  equation  (6.11)  into  the  form 

,  iafi2  Ve-»>  .  . 

3irrs  -  T^r  ’  S^SST !  » *  !■ 2-  3  «?-n 

*  n1 

The  right  hand  side  of  (9.1)  can  be  evaluated  from  the  power  spectral 
density  measurements  (notice  the  necessity  for  having  absolute  calibration 
of  electrie  and  magnetic  field  values)  and  from  the  previously  estimated 
8  and  0.  Since  the  solution  of  the  three  equations  (9.1)  results  in 
three  sets  of  four  possible  angles,  the  corresponding  angles  from  each 
set  were  averaged.  Table  I  shows  estimated  distance,  0,  and  possible 
azimuths  for  &  =  0°,  5°,  10°,  and  20°.  Al?o  displayed  in  the  table 
are  the  minimin  values  of  the  ft  function  (see  equation  6,15)  correspond¬ 
ing  to  the  "best"  values  of  e  for  a  given  6. 
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Table  I:  Source  location  estimates  based  upon  the  spectra  of  Figure  (9-1). 
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The  coil  used  for  the  magnetic  field  measurements  near  Kingston  is 
oriented  with  its  axis  in  the  geomagnetic  east-west  direction  and  the  0 


values  in  Table  1  were  therefore  corrected  for  the  magnetic  deviation  in 
Rhode  Island  (14°  west).  All  G  and  $  values  in  Table  I  are  close  to  the 
known  location  of  the  source  if  it  is  assumed  that  the  closest  azimuth 
estimate  is  the  correct  one.  The  maximun  distance  error  is  only  5° 

(about  300  miles)  and  the  maximum  bearing  error  is  7°.  The  fact  that  the 
0  and  0  values  which  are  closest  to  the  known  source  position  yield  the 
largest  Q.  parameter  does  not  seem  to  be  consistent  with  the  d«velopcien. 
in  chapter  VI.  In  chapter  VIII,  however,  it  was  shown  that  additive 
noise  can  cause  small  errors  in  the  location  estimates  and  this  is 
probably  the  case  here.  The  value  of  A  which  appears  to  provide  the  best 
agreement  between  theoretical  and  experimental  ratios,  as  indicated  by 
minimum  Q,  is  A  *.  20“  which  is  unrealistically  large.  Figures  (?.l) 
and  (7*2)  show,  however,  than  for  sources  located  at  distances  near  30°, 
the  spectral  ratios  are  nearly  independent  of  A  so  that  the  ratio  method 
can  provide  very  little  information  about  the  areal  extent  of  the  source. 
The  preceding  example  shows  that,  at  least  at  this  particular  time  (June  6, 
19'->7)«  the  background  noise  estimates  suggested  in  chapter  8.2  are  much 
too  large  and  correct  location  of  a  major  thunderstorm  region  is  clearly 
possible. 

9.2  Location  of  Two  Simultaneously  Active  Sources 

The  power  spectra  of  the  vertical  electric  field,  the  geo¬ 
magnetic  north-south  magnetic  fic-ld  and  the  geomagnetic  east-west 
magnetic  field  measured  in  central  Rhode  Island  on  January  1,  1970 
between  1720  and  1729  GMT  are  shown  in  Figure  (9.2).  Employing  the 
methods  developed  in  section  (7.2),  the  distances  from  the  receiver  to 
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each  of  the  two  sources  were  estimated  tc  be  0^  «  120°,  0p  =  lb2°  and 
the  relative  strength  parameter,  6,  was  estimated  to  be  4.22.  The  pos¬ 
sible  aslmuthal  locations  of  each  source  were  determined  from  the  ratios 
n  *  1,  2,  3  (see  section  7.3)  •  The  eight  locations  listed 
below  are  derived  from  these  calculations. 

Source  A 

(1)  13°N,  78°E  Southeastern  India  (2300  local  time) 

(2)  43°S,  20°E  South  African  Coast  (1900  local  time) 

(3)  64°S,  l4o°W  Southern  Pacific  (0700  local  time) 

(4)  2°N,  1>2  New  Guinea  (0200  local  time) 

Source  B 

(1)  28°S,  102°E  West  Australian  Coast  (2400  local  time) 

(2)  53CS»  90°E  Southern  Indian  Ocean  (2300  local  time) 

(3)  60°S,  120°E  Southeastern  Indian  Ocean  (0100  local  time) 

(4)  30°S,  125°E  Southwestern  Australia  (0200  local  time) 

Location  (3)  for  Source  A  can  be  eliminated  since  it  is  too  far 
south  (only  about  200  miles  from  the  Antarctic  Circle).  Locations  (2) 
and  (3)  for  fouree  B  are  also  unlikely  locations  for  thunderstorm  loca¬ 
tions.  It  can  be  concluded  that  Source  A  is  either  in  Southeastern 
India,  off  the  South  African  Coast,  or  in  New  Guinea;  that  Source  B 
lies  in  Western  Australia,  and  that  the  electrical  intensity  of  Source  B 
(total  power  output  of  the  thunderstorm  region)  is  4.22  that  of  region  A. 

It  is  generally  believed  that  the  worlds  major  thunderstorm  regions 
are  concentrated  in  tropical  and  subtropical  land  masses.  The  above 


results  do  not  contradict  this  statesent  and  curves  by  itaith  (1961)  giv¬ 
ing  the  average  frequency  of  occurrence  of  thunderstonas  throughout  the 
world  for  "Dec esrber-J anuary-February"  (Smith,  fig.  9*30)  indicate  rela¬ 
tively  high  thunderstora  activity  for  Southeastern  India,  the  South-Afri- 
can  -oast,  New  Guinea  and  Western  Australia.  Specific  surface  observations 
for  January  1,  1970  are  not  yet  available  to  us  at  the  present  time. 

As  a  further  example,  the  analysis  of  the  data  fcr  1333-1343  OT  July  29, 
1969  resulted  in  the  following  source  locations  (excluding  four  points  at  ex¬ 
treme  southern  or  northern  latitudes,  55°S,  T5°S,  55°N* 

Source  A 

(1)  28°N,  55°W  West  Central  Atlantic  (0900  local  time)  /  _  -o» 

(2)  25°N,  78°W  Bahama  Islands  (0900  local  time)  1  0  J 

Source  B 

(1)  1°S,  96°E  Sumatra  (2000  local  time)  ,  _o. 

(2)  5°S,  l40°E  New  Guinea  (2200  local  time)  (  5  ' 

Although  southern  hemisphere  surface  observations  were  not  available 
for  this  period  information  was  available  from  the  U.S.  Weather  Bureau 
"Northern  Hemisphere  Surface  Charts."  At  1200  UT— one  and  one  half  hours 
before  the  ELF  obseivations — the  weather  conditions  in  the  estimated  storm 
region  (20°-7O°N,  73°-83°W)  included  cumulonimbus  and  towering  cimulus 
clouds,  roughly  half  overcast,  and  wind  out  of  the  southeast.  At  1800  OT 
the  amount  of  cumulonimbus  clouds  had  greatly  increased,  precipitation  was 
oc curing  and  thunderstorms  were  observed  in  the  northwest  corner  of  the 
source  region  (northern  Florida).  It  is  thus  likely  that  lightning  ac¬ 
tivity  developed  soon  after  1200  OT  and  then  moved  to  the  northwest  of  the 
"source  region"  under  the  influence  of  the  prevailing  winds.  It  is  also 
possible,  however,  that  lightning  activity  was  already  present  at  1200  OT— 

although  not  indicated  on  the  "surface  chart"— since  only  seven  observa- 

6  2 

tion  stations  in  the  source  region  of  10°  kin  area  reported  weather 


phenomena  on  July  29  at  1200  UT. 
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X.  SUGGESTIONS  POR  FUTURE  WORK 

10.1  Possible  Improvements  in  the  Thunderstorm  Location  Method 

There  are  two  limitations  to  the  amount  of  information  about  thunder¬ 
storm  locations  and  characteristics  which  can  be  obtained  from  a  single 
receiving  site.  These  are:  (l)  the  ambiguities  in  tne  estimates  of 
source  azimuth  which  were  explained  in  chapter  VI  and  (2)  the  effects 
of  equatorial  background  noise  which  could  possibly  (see  chapter  VIII) 
be  large  enough  to  produce  serious  errors  in  the  so’irce  location  esti¬ 
mates. 

The  effects  of  both  of  these  difficulties  could  be  reduced  by  the 
utilization  of  simultaneous  measurements  tak  n  at  several  appropriately 
positioned  receiving  stations.  It  has  been  shewn,  for  example,  in 
chapter  VI  that,  if  the  background  noise  level  is  sufficiently  low,  data 
from  a  single  receiving  station  are  sufficient  to  determine  the  source 
distance;  bearing  information  is,  however,  ambiguous  since  choice  among 
four  possible  angles  is  arbitrary.  On  the  other  hand,  if  spectra  from 
two  widely  separated  stations  are  known,  only  two  possible  source  inn?.  - 
tions  exist  at  the  points  of  intersection  of  the  constant  distance  loci 
from  the  two  receivers.  Access  to  measurements  from  three  stations  would 
remove  all  ambiguity. 

The  preceding  considerations  assume  that  data  from  each  station  are 
processed  individually  to  obtain  possible  source  locations  and  that  the 
final  estimates  are  compared.  This  is  analogous  to  separately 


pressing  the  sign* Is  from  each  element  of  an  antenna  array j  we  are 
clearly  losing  Information  by  doing  this.  A  better  method  would  be  to 
collectively  analyze  the  raw  multi-station  data  using  techniques  which 
would  be  extensions  of  the  source  location  .methods  described  in  this 
report.  It  would  seem  likely  that  reliable,  continuous  monitoring  of 
global  thunderstorm  activity  could  be  obtained  in  this  uay. 

The  use  of  multiple  receiving  sites  could  also  improve  the  noise 
rejection  properties  of  any  thunderstorm  location  method.  If  the  stations 
were  placed  so  that  the  background  noise  spectra  were  different  at  each 
receiver  (because  of  different  geographical  locations  relative  tc  the  noise 
source) ,  the  effect  of  noise  contamination  could  probably  be  reduced  by 
cross -correlation  techniques.  It  might  also  be  possible  to  estimate  the 
noise  spectrun  and  analytically  remove  it  from  the  received  spectra.  In 
short,  there  is  no  doubt  that  multiple  station  measurements  could  pro¬ 
vide  more  reliable  location  estimates. 

10.2  Estimation  of  the  Ionospheric  Conductivity  Profile 

As  was  pointed  out  in  Chapters  IV  and  V,  Schumann  spectra  depend 
up .n  both  the  location  and  characteristics  of  the  source  and  the  prop¬ 
erties  of  the  cavity  boundaries,  particularly  the  ionosphere.  It  was 
shown  that  the  relative  amplitudes  of  the  peaks  in  the  resonant  spectrum 
provide  a  measure  of  source  location  which  is  nearly  independent  of  the 
properties  of  the  ionosphere  (see  Figures  5.8  to  5.1^).  On  the  other 
hand,  the  exact  spectral  position  of  the  resonant  frequencies  depends 
about  equally  on  the  source  location  and  the  ionospheric  conductivity 
profile  (see  Figures  5.15  to  5.20). 


It  is  possible  t.ier_fore  to  first  determine  the  location  (0)  and 
extent  (A)  of  the  source  from  ratios  of  spectral  peak  powers  and  then 
to  employ  the  resonant  frequency  values  to  draw  conclusions  concerning 
th»  conductivity  profile  •  f  the  ionosphere.  The  Schumann  spectra  are 
affected  by  regions  of  ionization  extending  as  low  as  JO  to  40  km  above 
the  earth  surface  (Jones,  1967)  and  therefore  should  provide  information 
which  is  difficult  to  obtain  by  other  means.  For  example,  if  it  is  de¬ 
termined  from  the  ratios  that  0  =  100°,  A  =  20°  and  if  the  first  reson¬ 
ant  frequency  in  the  electric  field  is  at  8  Hertz,  then  Figure  (5.18) 
suggests  that  the  "Pie rc e-Cole"  profile  is  a  more  appropriate  ionosphere 
model  than  the  Galejs  profile. 

A  refinement  of  the  methods  for  both,  source  location  and  deter¬ 
mination  of  the  Ionosphere  profile  would  involve  an  iteration  procedure. 
This  would  require  the  determination  (from  the  observed,  experimental 
ratios)  of  the  location  and  extent  of  the  source  assuming  firot  a  reason¬ 
able  average  ionosphere  profile  such  as  the  "Pierc e-Cole"  profile  shown 
in  Figure  (5.7) •  Then  the  exact  values  of  the  observed  resonant  fre¬ 
quencies  and  the  estimated  source  positions  would  be  used  to  obtain  a 
"better"  ionosphere  model.  This  model  would  then  be  used  again  to  obtain 
a  better  estimate  of  the  source  location.  The  procedure  would  be  continued 
until  further  iterations  produce  no  change  in  the  ionosphere  profile. 

10.3  Determination  of  Relative  Contribution  of  Major  Thunderstorm 
Regions  to  World  Electrical  Activity 

As  indicated  at  the  end  of  chapter  VIII,  verification— or  lack  of 
verification — of  ELF  thunderstorm  locations  by  ground  or  satellite 
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observations  should  provide  information  about  the  relative  importance 
of  the  major  thunderstorm  regions  as  lightning  generators  in  comparison 
with  the  totality  of  widely  distributed,  small  thunderstorms.  A  sys¬ 
tematic  program  of  comparing  the  results  of  Elf  spectral  analysis  with 
surface  observations  would  therefore  provide  information  about  the 
diurnal  and  seasonal  variability  of  the  ratio  of  these  two  inputs  to 
world  electrical  activity. 

10.4  Investigation  of  Localized  Ionospheric  Perturbations  on  Schumann 
Resonances 

It  has  been  indicated  in  chapter  VHI  that  the  Schumann  spectra 
are  sometimes  contaminated  by  l/f**  noise.  One  possible  source  of  such 
noise  may  be  hydromagnetic  wave  generation  processes  in  the  Magnetosphere 
or  the  Solar  Wind.  The  l/f®  noise  could,  however,  be  also  a  manifestation 
of  near-field  effects,  that  is  superposition  of  TE  and  higher  order  TM 
modes,  produced  by  a  local  lateral  ionospheric  inhomogeneity,  or  it  could 
be  due  to  random  current  fluctuations  in  the  local  Ionosphere. 

To  determine  the  source  of  l/f®  noise  comparison  of  Elf  measurements 
with  Ionospheric  and  other  geophysical  data  such  as  cosmic  ray  count, 

10.7  cm  flux  and  micropulsation  activity  is  necessary.  Detailed  investi¬ 
gation  of  analytical  models  describing  the  effect  of  local  Ionospheric 
perturbations  would  also  be  required.  Madden-Thompson  (1965)  have  pointed 
out  that  the  global  effect  at  ELF  of  local  Ionospheric  perturbations  would 
be  small;  higher  order  modes  may  nevertheless  be  appreciable  at  ELF  within 
several  hundred  kn  of  an  inhomogeneity  or  discontinuity.  For  this 
reason  it  would  be  particularly  useful  to  study  the  ELF  spectra  corres- 


ponding  to  known  Ionospheric  discontinuities  such  as  ihcse  exis 
local  dawn  or  dusk  or  during  a  Solar  eclipse. 
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APPEtTuX  A 


Spectral  Bias  for  Rectangular  Data  Window 


The  quality  factor  for  resonant  cavities  is  defined  a3 
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where  U  is  the  energy  in  the  cavity  and  u)q  is  the  radian  resonant  fre¬ 
quency.  Equation  (A-l)  is  a  differential  equation  for  U  which  has  the 


solution 
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where  Uq  is  the  initial  3tored  energy.  Since  power  is  the  time 
derivative  of  the  energy,  the  power  behaves  in  a  similar  manner. 
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From  equation  (A -3)  the  fields  themselves  oscillate  as 

-  t  ioj  t 
E(t)  =  E  e  ° 


In  the  frequency  domain,  the  behavior  is  described  by 
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li  there  is  no  energy  in  the  cavity  for  t  <  0.  Integrating  (A-5), 


Vjui. 


in  the  region  of  each  resonance  and  determine  the  importance  of  the 
additive  bias  term  S"(f).  The  power  and  the  power  spectrum  are 

approximately  related  by 


P(f)  =  S(f)Af 


(A-10) 


so  that  the  n.rmalized  additive  bias  term  for  spectral  density  esti¬ 
mates  is 


The  frequency  variation  of  this  bias  tern  in  the  region  about 
each  of  the  first  three  cavity  resonances  is  shown  in  Figure  (A-l)  for 
a  frequency  resolution  f  =  1  Hz.  It  can  be  seen  from  Figure  (A-l)  that 
near  the  first  resonance  the  additive  bias  of  the  estimates  is  about 
-l8  percent  of  power  spectral  density  at  that  point.  This  rather  large 
error  can  be  reduced  by  decreasing  the  resolution  to  0.5  Hz  or  0.25  Hz, 
which  reduces  the  Dias  to  about  -4  percent  and  -1  percent  respectively. 
From  the  preeeedirg  discussion  it  is  clear  that  if  resolutions  of  0.5 
or  less  are  used,  ;he  bias  produced  by  employing  a  rectangular  data 
window  is  negligible. 

An  intuitive  argument,  verifying  quali ;atively  the  mathematical 
results  presented  above,  is  provided  by  Figure  (A-2)  where  a  very  high 
Q,  resonance  is  shown.  The  procedure  outlined  in  Chapter  II  uses  the 
average  power  over  a  bandwidth  Af  as  the  estimate  of  the  power  spectrum 
at  the  center  of  the  band.  In  Figure  (A-2)  the  average  power  in  the 
bandwidth  Af  is  schematically  represented  by  point  A;  m  the  other  hand, 


the  actual  power  at  the  center  of  the  band  is  given  by  point  B.  It  is 
clear  that  approximating  the  power  apectrun  by  the  average  power  in  a 
band  Af  can  be  very  misleading  if  the  spectrum  is  rapidly  varying  with 
frequency. 

This  problem  can  be  eliminated  by  decreasing  Af ,  but  this  will  in¬ 
crease  the  variance  of  the  power  spectral  density  estimates  as  is  indi¬ 
cated  by  equation  (2.19).  There  is  then  an  optimum  range  of  Af  which 
must  be  used!  Af  must  be  chosen  large  to  minimize  the  variance  of  the 
estimates;  but  it  should  be  made  small  to  minimize  the  bias  of  the 
estimates.  This  appendix,  together  with  Chapter  II  has  shown  that  a 
Af  of  0.25  Hz  or  0.5  Hz  provides  acceptaole  levels  of  estimate  bias  and 
variance. 

It  is  not  to  be  construed  from  the  above  that  the  lower  the  Q,  of 
the  resonances  (or  equivalently,  the  less  the  frequency  variation  of 
the  spectrum)  the  better  the  information  contained  in  the  spectrum.  Ob¬ 
viously  if  the  spectrum  is  white  (ft  =  0)  the  bias  is  zero  but  there  is 
no  information  provided  by  the  spectrum.  The  purpose  of  this  appendix 
is  merely  to  demonstrate  that  for  a  frequency  spectrum  of  the  Schumann 
type,  where  the  quality  factors  of  the  resonances  are  of  the  order  of 
4  or  less,  a  frequency  resolution  of  0.25  Hz  or  0.  5  Hz  gives  good  esti¬ 
mates  of  the  actual  spectrum. 


APPENDIX  F* 


Linear  Stationary  Systems  with  Uncorrelated  Inputs 

If  a  linear  system  is  driven  by  a  source  at  radian  frequency  u, 
the  Fourier  transforms  of  the  output  y(t)  and  input  x(t)  are  related  by 
the  transfer  function  H(jw)  as 

Y(Ju>)  =  H(ju)  X(jw)  (B-l) 

H(  jui)  is  the  Fourier  transform  of  the  impulse  response  of  the  system. 

It  can  be  shown  (Davenport  and  Root,  1938)  that  the  output  power  spec¬ 
trum  is  given  by 

Sy(u»)  =  |H(»j2  Sx(w)  (B-2) 

If  we  generalize  this  to  a  multiple  input  system,  the  single  output 
y(t)  is  related  to  the  N  inputs  X^t);  k  *  1,  2,  . ..N  by 

N 

y(t)  =  J'  y  (t)  (b-3) 

n=l 

where  yk(t)  is  the 

system  and  the  other  inputs  are  zero.  The  y^(t)  are  related  to  the  x^t) 

by  the  convolution  integral 

+00 

yk(t)  =  j  hk(T)xk(t-T)dT 


response  when  only  the  input  x^t)  is  driving  the 


(B-4) 
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or  equivalently 


\(jw)  -  H^Jw)  yjw) 


(B-5) 


where  n^(Jw)  and  h^(t)  are  the  transfer  functions,  in  the  frequency  and 
time  domains  respectively,  between  the  kth  input  and  the  output  when 
all  other  inputs  are  zero. 

If  the  inputs  x^,  ...  are  random  signals,  the  output  power 

spectral  density  is  given  by  (Davenport  and  Root,  1938) 


N  N 

VW)  -  l  ZHn*(jw)Hk(j“)Snk(w) 

k>l  n-1 


(B-6) 


whore  3nlc(ai)  is  the  cross  spectral  density  between  the  inputs  xQ(t)  and 


x^t).  SjjjjCw)  be  expressed  as 
+“ 

**W-J  A  x, 


(B-7) 


where 


R  (t)  is  the  cross  correlation  function  of  x  (t)  and  x.  (t), 

X  X.  II  it 

n  k 


If  the  inputs  are  uncorrelated,  8^(01)  *  0  and 


s  (f)  -  Z  IVJ“)l2sx(“) 


(B-8) 


which  the  reader  should  compare  with  eq.  (B-2). 

A  statement  of  the  result  given  by  (B-8)  is  as  follows  1  The  power 
spectral  density  of  the  output  of  a  multiple  input  linear  system  driven 
by  uncorrelated  inputs  is  equal  to  the  sum  of  the  power  spectral  den¬ 
sities  due  to  each  input  considering  all  other  inputs  to  be  zero. 
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He  can  apply  the  results  of  this  appendix  to  the  problem  of  light¬ 


ning  excitation  of  the  earth-ionosphere  cavity  by  aaswlnf  that  the 
signals  frcai  the  thunderstorms  are  stochastic  in  nature  and  that  the 
different  lightning  strokes  are  Independent.  Hence,  the  signals  from 
each  discharge  are  uncorrelated.  Ra enter  (1961)  has  pointed  out  that 
this  approach  Is  mathematically  equivalent  to  the  analysis  of  "shot 
noise"  in  vacuum  tubes  where  it  is  usually  assumed  that  the  probability 
of  the  mission  cf  an  electron  from  the  cathode  is  independent  of  the 
number  of  electrons  emitted  in  the  past  (Davenport  and  Root,  1958) . 
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Figure  5.2.  ELF  phase  velocities, 
D.  L,  Jones  (1967). 


a)  Piero*  and  Cole  ionosphere  profile  emitting  data  below  5 6  km 

b)  Desks  ionosphere  profile  (B.  R.  May,  J.  Ataoaph.  Terr.  Phya.,  28. 

55Ji  1966) 

c)  Cole  and  Pierce  (1965)  ionosphere  profile  including  heavy  ion  region 
down  to  350  km  altitude 

d)  Profile  (b)  continued  to  lower  altitudes  by  profile  (c) 

e)  Chapman-Jones  (1964)  two-layer  profile. 
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Figure  5.3.  ELF  attenuation, 

D.  L.  Jones  (1967). 


a)  Pierce  and  Cole  ionosphere  profile  omitting  data  below  56  km 

b)  Deeks  ionosphere  profile  (B.  R.  May,  J.  Atmosph.  Terr.  Phy i . ,  28, 

553,  1966) 

c)  Cole  and  Pierce  (1965)  ionosphere  profile  including  heavy  ion  regi 
down  to  30  km  altitude 

d)  Profile  (b)  continued  to  lower  altitudes  by  profile  (c) 

e)  Chapman-Jones  (1964)  two-layer  profile. 


Figure  Ratio  of  resonant  peaks;  sources  at  9  -  L  to  9  +  A. 

Results  of  approximation,  equation  (5. 19) ,  and  results  of  numerical  Integration 
employing  equations  (5-26)  and  using  v  based  on  Cole-Pierce  profile. 


Ratio  of  resonant  peaks;  magnetic  field;  sources  at  8  -  A  to  8  +  a 
Results  of  approximation,  equation  (5.20)  and  results  of  numerical 
integration  of  (5.2?) ;  v  based  on  Cole-Pierce  profile. 
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Figure  5.16.  Variation  of  magnetic  field  second  resonant  frequenoy  with  source  distance 
for  two  different  ionosphere  profiles. 
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5.20.  Variation  of  electric  field  third  resonant  frequency  with  source  distance 
for  two  different  ionosphere  profiles. 
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Figure  6.1.  Flowchart  describing  method  for  determining  best  distance 
estimate  for  two  measured  ratio  values. 


figure  8.1.  East-We3t  magnetic  field  spectrum  in  Rhode  Island  at  0400  EST  August  30,  1967. 
An  illustration  of  cosmic  noise  contamination. 
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Figure  R.4.  Estimated  electric  field  spectra  for  uniform  lightning  distribution: 

summer  lCf'S  to  60°N,  winter  40°S  to  10°N.  fall-spring  30° S  to  jW°N. 


Figure  8.7.  E-W  magnetic  field,  March  2j5,  19$7 


Figure  P.°.  Vertical  electric  field,  August  14,  19*59*  1400  EST 
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